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ABSTRACT 


An experimental apparatus to measure the particle density 
distribution in the freeboard of an atmospheric fluidized bed was 
designed and constructed. The density versus height measured by 
the sampling apparatus gives a similar exponential decrease as 
previous investigations have found. 


A particle trajectory model is developed which calculates the 
height and particle density distributions above the bed surface of 
an atmospheric fluidized bed. The parameters input to the model 
are the superficial velocity, initial partical velocity, gas jet 
velocity and duration, and the particle size distribution of the 
bed mass. The model was evaluated using the experimental data for 
jet velocity, duration, and particle size. The predicted slope of 
the particle density versus height in the freeboard agrees with 
the experimentally measured slope within 20%. 


A sensitvity analysis using the trajectory model, resulted in 
a determination of the particle distributions in the freeboard of 
a fluidized bed as affected by varying the input parameters to the 
trajectory model. The most significant effects were achieved when 
the jet velocity or duration was altered. 
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CHAPTER.. I 


INTRODUCTION 


Fiuidized beds have been used in industry for many years. 
They have been used to mix and dry particulate materials and are 
the princigle process in catalytic cracking plants. In the past 
decade or so, the use of fluidized bed combustors for power 
Qeneration has become a source of major interest. Prototype coal 


burning beds have already been built which are comparable to 


یف سا 


existing coal plants. Fluidized bed combustors have the added 
benifit of low NOx, S05 and hydrocarbon emissions and the 
flexibility ot being able to burn a uide range of fuels ranging 


from refuse and high sulfur content coal to high grade fuels. 


A fluidized bed [Fig 1] is composed of a distributor through 
uhich an air flou is introduced through thousands of small 
orfices. This air then passes through the dense zone of the bed 
which is comprised of amass of particles. The air E 
through the dense zone is maintained above the minimum 
fluidization velocity (Umf) during normal operations. At 
veloci*ies ecual to or greater than Umf, the frictional force 
(Drag) of the air flouing pest a particle is equal to the weight 
of the particle. Under these conditions, the particle mass 
behaves very much like a fluid. It will maintain a horizontal 


surface if the container is tilted, flow out of holes in the 
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container, and has a pressure drop across any section of the bed 
apprxamately equal to the weight of air and particles in the 
section [1]. At velocities above Umf, the excess air will 
coalesce and form small voids or bubbles as it procedes towards 


the surface of the fluidized bed. 


As bubbles break at the surface of the bed, the solid 
particles are thrown up above the bed surface and are entrained by 
the upward flowing gas stream. This zone above the bed surface is 
the freeboard zone. In the freeboard, some particles are carried 


by the gas flow far atove the bed surface and are removed from the 


— 
ኤኤ... 


fluidized bed (elutriated) , while the remainder fall back to the 
bed. In general the amount of bed solids suspended in the 
freeboard (entrainment) decreases exponentially along the freeboard 
height. This distribution is similar to that of the 
Maxwell-Boltzmann distribution for the case of discrete energy 


states as it applies to the Law of Atmospheres [2]. 


N(z) = No expt-mgz/kT} (1) 


Research in the area of entrainment by Lewis et al. [3], Zenz 
and Weil [4] and others has resulted in the following correlation 
for entrainment as a function of gas velocity and freeboard height 


for small-particle beds [1]. 
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2 
= Bexpt - ((b/Uo) + aH J} 
At Uo (Z) 





where: 
F = Entrainment 

At = Area of bed 

Uo = Superficial velocity through bed 
= Particle dependent constant 
Particle dependent constant 
= Particle dependent constant 
= Height of freeboard 


rag 0 0 
ዘ 


At some point above the bed surface, the quantity of 
entrained particles becomes constant. At this point, the free 


fail velocity of the remaining particles is equal to or less than 


ሦ.. 
~ 


the uniform superficial operating velocity. The height at which 
the entrainment becomes constant is called the transport 


disengagement height (TDH) (11. 


The particles that are thrown above the bed are affected in 
the freeboard region by hydrodynamic parameters such as: bubble 
size, bubble frequency, fluidizing velocity, height above the bed, 
particle size, particle density, column diameter [51, and 
Daffles [8]. The intermittent high velocity bursts of gas which 
occurs when a bubble bursts, imposes a flucuating and highly 
irregular time dependent velocity profile over the cross section 
of the bed surface. At successively higher levels above the bed 
surface, this velocity profile becomes more and more uniform until 
at the TDH, the flou is at the uniform superficial operating 


velocity (Uo) [4]. 
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Until recently, little attention was paid to the 
understanding of the freeboard reactions for large particle beds. 
However, Gue to recent work in fluidized bed combustors, the 
extent of particle and fines loading in the freeboard has been 
shown to significantly affect the 5ዑ7 absorbtion, NOx reduction, CO 
emission. In general, the last S to 10% of the E te will 
burn in the freeboard. It uas shoun that the fine sorbent 
particles entrained into the freeboard will enhance sulphur 
capture and that entrained char particles will react with NOx and 


reduce its emission [7]. Fines reiniection has been shown to 


significantly. increases the fine particle concentration in the bed 


Dn 


and in the freeboard with the consequence of further enhancing 
char oxidation. However, this can result in overheating in the 
freeboard region and excess 502 and NOx emission [7]. The end 
result is that the potential for fluidized bed combustion power 
plants is enhanced by their ability to burn high sulphur content 
fuels and maintain lou 30, emissions. Further research in the 
area directly above the bed surface is required to properly model 


the reactions occuring within the freeboard. 


Extensive studies on entrainment rate and elutriation have 
been made with numerous correlations, some of which are proposed 
in (1, 3, 4, S, 6, 7, 8). However, most of the reported work on 
entrainment from fluidized beds has been carried out with either 
a closely sized fraction of particles or a mixture of two such 
fractions. Virtually all of this work has been conducted on bench 


scale or catalytic cracking fluidized beds. The results of these 
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studies have been shown to produce much lower entrainment rates 
than full sized beds or tend to operate in the slugging condition. 
Therefore, the entrainment rates and transport disengagement 
heights (TDH) for fluidized beds are generally estimated from 
empirical or semiempirical correlations obtained from this data 
and most of them show extreme discrepancies between different 
experimental results. -Extrapolation of these empirical 
correlations usually leads to strange results [5] with 


discrepancies which can vary by two orders of magnitude. 


— The lack of good correlations stems mainly from the 


— 
— 


difficulty in obtaining accurate entrainment rate data. Most of 
the data is based on pressure measurements at incremental heights 
in fluidized cracking plant type beds [111. The effect of wall 
loading by particles and the actual relation — pressure and 
particle concentration is considered to be major sources of error 
when using this method with large particles. As a result, none of 


the correlations are widely accepted as giving accurate 


predictions [8]. 


A complete model of the entrainment process from fluidized 
beds must take into account all the mechanisims involved within 
the process. The arrival of bubbles at the bed surface, ejection 
Of particles from the dense-phasec bed into the freeboard region 
as the bubbles eruot, particle-particle interactions, and the 
trajectories of ejected particles are all important [8]. Much 


work has been done concerning bubble grouth, velocity, volume, 


3 
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etc. and their behavior is fairly well understood. 


The mechanism of solids ejection at the bubbling bed surface 
is still not well understood. The origin of ejected particles is 
reported to be primarily due to two sources. The particles which 
have been lifted by the bubble wake and thrown upwards following 
the bubbie burst at-the surface-is the-first source. This theory 
is supported by work done by George and Grace (8] who performed 
experiments which concluded that the vast majority of the ejected 


particles did not originate from the surface layers but from 


‚Bubble wake pick up.-..Work done by Page and Harrison [6] also 


— 
— 


appears to agree with this, The second theory suggests that the 
ejected particles originate at the nose of the bursting bubbles 
and are thrown outward when the bubble breaks. Research by Rowe 
and Partridge [8] and Glicksman et al [12] have 00 this second 
mechanisim as being the dominate particle ejection source and thus 
supporting this second theory. Their work has also shown: that 
under the conditions in which 2 bubbles coalesce just below the 
surface of the bed, the jet of gas produced can result ina 
significant amount of particles being ejected from the wake of 


the first bubble. 


The effect of multiparticle interactions have been for the 
most pert ignored except by Peters and Prybylouski [13]. The 
motion of any individual particle is influenced by the presence 
of other particles, i.e., through direct particle-particle 


interactions and deviations in the fluid drag force. The major 
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drauback of their work is that the paper compares their theory 


with only a single set of experimental results [3]. 


Studies to model the trajectories of particles in the 
freeboard have been conducted several times. The work of Walsh 


et al_{7], George and Grace [8], and Peters and Prybylouski [13] 


- .name just-a feu of the latest-efforts. All of these studies 


relied upon experimental data to develope their theories. 
However, to check the accuracy of their theories, more 


experimental data is required. 


As of yet, none of the entrainment models available can be 
incorporated into fluidized bed combustion models with sufficient 
accuracy to warrent their use. This is due to a lack of 
experimental information on entrainment rate as = function of the 
complete fluidization parameters of the bed to test the models 
with. As a result, the purpose of this study has been to obtain 
particle density distributions above a cold atmospheric fluidized 


bed containing a continuous particle size distribution [Appendix 
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CHAP FER I Y 


PARTICLE SAMPLING APPARATUS 





Design Alternatives 


There are many methods available for determining particle 


distributions in fluid flows. The more commonly used methods are: 


1) Catching mechanisms 
۲: سوستج‎ 2) Trapping mechenisms 
3) Radiation attenuation measurements 
4) Optical measurements 


5) Capacitance and Inductance measurements 


Catching mechanisms are passive devices. That is, particles 
are captured merely by the presence of the catching mechanism in 
the fluid flow containg the particles to be sampled. The data 
obtained using this method is position dependent and produces 
average values for the particle flux loadings. These catching 
mechanisims are also limited in that they can only catch particles 
with particle fluxes traveling in a single direction. The device 
used by Walsh et al [12] only captured falling particles while the 
device used by George end Grace (8] required the uouard moving 
particles to deflect off of a baffle surface and fall into a 


collecting trough. 


IN 
IN 





Trapping mechanisims, unlike catching mechanisms, are active 
particle samplers. Their operation involves the trapping and 
isolation of a finite volume of the fluid flow at a specific 
period in time. This sampling technique produces time dependent 
as well as position dependent data. This will allow correlations 
between bubble eruption and particle density to be made using 
muitipie bubble conditions cather that single bubble capture. As 
the number of random samples taken by this method increases, the 
average value of this data will approach that of the catching 
mechanism. Trapping mechanisims also capture particle fluxes 
traveling in multiple directions.= This ability reduces the error 
inherent in measuring oniy the downward or only the upward 
particle flux. The apparatus used in this paper is a trapping 


mechanism. 


Attenuation of nuclear particles from a radioactve source can 
be used to give average particle density distributions across a 
suspension. However, this method is not adaptable to density 
determinations at a point. An average time dependent density 
determination can be achieved with this method. Another drau-back 
of this method is the radiation hazards involved with the use of 


nuclear particles. 


Optical density determinations consists of two seperate 
methods. The first method uses a very small light beam which is 


eclipsed by the transition of a particle through it. A related 
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method uses the absorbtion and scattering of a someuhat larger 
light beam to correlate the change in light intensity with 
particle density. This method has been used frequently in the 
study of aerosols but reqwires complicated and intricate equipment 
(141. The second method involves high speed photographs of a 
small volume of space. This method cannot be used when the 

eur. -Particle .density-is.so.large that multiple particles eclipse each 
other frequently enough to produce unacceptable error. This is 
the case when the probe height above the bed is less than 7-15 cn 


Eso 1n). 


--- Measurements at a point can also be made by inserting either 
a toroidial inductor or a parallel plate capacitor in the flow. 
The presence of the particles changes the permeability and thus 
the inductance of the inductor, or the dielectric strength and 
thus the capacitance of the capacitor. The drau-backs of these 
methods involves the unknown effects of particle velocity and 
external particles on the inductor and charge transfer to 


particles from the capacitor [14]. 


Apparatus Requirements 


The goal of this study was to determine the density 
distribution of particles above an atmospheric fluidized bed with 
particle velocities of up to 12 meters per second. The particie 
Size distribution of material ejected from bubbles is required for 


particle trajectory calculations. A correlation between the 
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average density and the density present immediately after a bubble 
bursts from the bed surface was also of interest. These 
requirements dictated that the method used for measuring densities 


have the following capabilities: 


1) Measure densities with good spatial resolution. 


2) Measure censities at specific moments in time. 


- 


— 


-:-. 3) Obtain particle size information. 
4) Operate under extremely dirty conditions. 


5) Easy sample removal from bed. 


6) Remote operation of sampler. 


The radiation attenuation and inductance/capacitance methods 
can not determine particle sizes. Therefore, these methods uere 
no longer considered as possible measurement alternatives. 

Because the optical methods are either not reliable at small 
heights above the bed or their use is too complex, they uere not 
used. Catching devices, although simple to use, do not have the 
ability to measure data at specific points in time and determine 
particle density loading in space. As a result, the determination 
to use a trapping mechanism as the method of measuring particle 


densities was made. 
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Apparatus Design 


General Design Criteria 


The following criteria was used in determining the design of 


the trapping device. 


1) 


2) 


4) 


The closure time of the trap was chosen to be equal 
to the time for a particle with a velocity of 18 
m/s to transit 1/10 the length of the sample 
container. This velocity is considered to be the 
upper limit of the particle velocity distribution 
present in the test bed, based on the uork of 


George and Grace [81. 


The apparatus must be capable of frequent sampling 
without requiring access to the sampling device 


itself. 


The samples trapped, must be easily accessable 
from outside the fluidized bed without interupting 


the bed conditions. 


The apparatus must be able to operate in the high 


particle flux environment of the fluidized bed. 





S) The actuation of the trapping device must be able 
to be accurately determined to allow correlation 


with other time resolved measurements. 


General Design 


The apparatus is shown in Figs. 2 and 3. A description and 
list of all components is given in Appendix H. The sample 
container is mounted on an extension arm to minimize the 


MMM‏ ہق 


—disturbance to the air flou around the sample trap caused by the 


- 
€ 


rest of the mechanism. The sample container is closed using tuo 
(2) paddle arms, one above and the other below. These paddle arms 
are attached to aluminum interfaces which are used to connect them 
to a rotary solenoid. The solenoid is used to suing the paddles 
over the sample trao and shut it. Not shown in these figures are 
the power supply for the solenoid, the vacuum system used to 
remove the particles from the sample trap and the water-proof 
nylon shell used to keep the particles from interfering with the 
operation of the solenoid. All of these systems are described in 


greater detail in the following sections. 


Sample Container 


To ensure that the sampling device had minimal effect on the 
fluid flow, the cross sectional area presented to the flow had to 


be minimized. This constraint required that the sample container 
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be seperated from the rest of the apparatus. This was benificial 
in the fimal design because it helped to reduce the apoaratus 


closure time. 


To ensure a good seal was achieved when the trap was shut, 
feit was used as a gasket between che sample container and the 
closure paddles:  Fig.:-4-shous.the final design used for the 
sampling mechanism. It is constructed of 1/16 inch aluminum with 


1/4 inch square stock used for the frame and mounting structure. 


Epoxy is used to seal the sides of the container. 


Eee Paddles 


After several iterations on paddle designs, it was determined 
that the paddle construction which offered the greatest stiffness 
for the least weight was a composite laminate. The paddle, shown 
in Fig. S, is made using a 0.4 in thick foam core with 1/32 inch 
thick Basswood laminations on the exterior. Hardwood (Maple) end 
pieces were used to provide a noncompressive connection between 
the foam paddles and the aluminum interfaces. The aluminum 
interfaces couple the solenoid shaft to the paddles. Epoxy was 


used to join the laminate materials. 


Actuator 


A rotating mechanism utilizing a rotary solenoid was chosen 


to shut the sample trap. A rotary solenoid was selected because 
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of the problems of maintaining low friction surfaces for sliding 
mechanisms in the presence of the particulate matter being 
sampled. This precluded the use of any sliding mechanism to shut 
the sample trap. The calculations found in Appendix A and 
Appendix C determined the size of the solenoid required to achieve 


the desired closures time. 


The 4S degree stroke solenoid was chosen to place the paddle 
arms far enough away from the sample trap, such that when 
de-energized it prevents interference with the particle flow. 


This stroke also minimized the area which must be clear of 


کے 
— 


z=. ——--obstructions-to-the travel-of the paddles. The solenoid operates 
at a 1/10 duty cycle power rating when initially actuated, 
providing the torque output shown in Fig. 6. After the solenoid 
has shut the sample trap, the solenoid is operated at a lower 
power rating, providing a holding torque of 5.S in-lbs. This 
decreased rating is necessary to prevent overheating of the 


solenoid. This assembly is encased within a nylon shell. 


Power Supply 


Fig. 7 is a schematic of the electrical system used to power 
the solenoid. Appendix H contains a list of all components used 
in the power supply. The power susoly plugs directly into a 
standard 11S volt AC line source. Switch Si is used to apply 
power to the solenoid Ml. The full wave bridge rectifier assembly 


converts the AC line voltage to OC. The rectifier assembly is 
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Fue. 8 Torque Output Curve for 
Rotary Solenoid. 
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protected by the arc suppresser Dl from the large voltage spike 
incuced in the solenoid coıl wnen the solenoid is actuated. 
Resistor R3 is used to reduce the current through the solenoid 
after the solenoid has been closed to prevent it from overheating. 
To initially shut the sample trap, full current is applied to the 
solenoid by shutting switch SZ, which bypasses resistor R3. 


Resistors RL and RZ form a voltage divider network to provide a 


lou voltage (5.5 volts AC) trigger source for an oscilloscope. 





Vacuum System 


To remove the particles in the sample trap, a vacuum system 
was developed. This system is shown in Fig. 8 with a list of the 
components given in Appendix H. A vacuum is produced by allowing 
air from a 10@ psi air source to flow through valve Vl into the 
venturi edcuctor Pl. The vacuum places a suction on the sample 
container CZ via a fine mesh screen. The purpose of this screen 
is to prevent particles from escaping the sample container. The 
suction is applied to the sample trap Cl through 1/4 inch polyflouw 


tubing. It is through this tube that the particles are removed 


— 
- 


ars the sample trap and collected in the sample container. An 
equalization and agitaion line is connected to the opposite side 
of the sample trap. This line serves two purposes. First, it 
ensures that the vacuum system does not pull les from 
outside of the sample trap. Second, it allows a flow of air to be 
introduced which stirs up the particles trapped inside. 5 
helos push them into the suction line and reduce the remaining 


pacticles to a minimum. 


Spoaratus Testing 


Two tests were run to determine the effectiveness of the 
system. The first test determined the closing time of the sample 
trap. The second evaluated the error from the loss of particles 


which were left in the sample trap by the vacuum system. 
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The procedure and method used to determine closure time for 
the sample trap is gıven in Appendix C. From the results of these 
tests given in Appendix C, the closure time was determined to have 
an average value of 1.44 milliseconds. This is equal to a 10 m/s 
particle traveling 1.44 cm or approximately 16 % of the sample 
trap length. The average particle velocity will be iess than 2 


m/s and will therefore introduce an average error of less than 3 1. 


Using the sample trap closure time, a dynamic analysis 
determined that the total time from initially applying power to 
the solenoid until it shut the sample trap is 42.6 milliseconds. 
مت‎ calculations are given in Appendix D. This actuation time 


is important for determining uhether a specific bubble's debris 


was within the vicinity of the trap at the time of closure. 


The procedure and results for the testing of the sample 
removal vacuum system are given in Appendix B. The average amount 
of particles lost by the vacuum system was determined to be 0.52 X 
of the initial sample placed in the trap. The maximum error was 
0.23 X. The error from the vacuum system is therefore considered 


to have an insignificant effect on the data obtained. 
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EOS TER IIT 


EXPERIMENTAL PROCEDURE 
Fluidized Bed Configuration 


[he M.I.T. atmospheric fluidized bed, in which the sampling 
device was used, is a model of the 29 MW atmospheric fluidized bed 
combustor prototype, jointly sponsored by the Tennessee Valley 
Authority and the Electric Power Research Institute. The 
tiuidized bed model is described in Lord et al [151 and Jones et 
al [151]. While using the sampling apparatus, a different heat 
exchanger Ms bundle configuration was used than is described in 
Jones et al [168]. The heat exchanger configuration used is shoun 


in Figs. 9 and 1@. 


The heat exchanger used during this work is made of 1.28 cm 
(0.5 in) O.D. tubing arranged in 4 rows of 22 pipes each. The 
tubes are aligned as shown in Fig. 10. Each pipe is spaced with 
a vertical center to center distance of 5.08 cm (Z in) and a 
horizontal center to center distance of 3.91 cm (1.5 in). The 
distance from the distributor to the center of the upper most tube 
is 27.52 cm (19.875 in). A 5.08 cm (2 in) spacing separates the 
front and back walls of the fluidized bed from the enc tubes of 
the bundel. The cross sectional area of the bed is 1.078 sq m 


(11.61 Sq Ft). The particulate material used in the bed is a 


40 





ገር). 





3 


80 


1.26 0 
° ی‎ Tubing 


Heat Exchanger Tube Design Shouing The 
Four (4) Rous of 22 Tubes. 
All Dimensions in cm. 


4] 


IA 


1 
مم 





ie 
Front ۳ | 
a لس‎ — 77.5 
_ نے‎ l 


Distributor 





6 


II - 
EN ۱ 
۱ د‎ 


= 9 سس 


Ui 
in 
6 
0 


د د ه د ه و ۵ ه٥‏ 5 د د د د د ۵ د ۵ 6 ۵ 6 6 د ه 


32999228 33293934. 0999233399953 ه‎ 
Side 


۱ م۸ 
a‏ 
ه ه ه ه ه ۵ ه ه ه ه ه ه و و ه ه ه ه ه ه ه ه ۵ ۱ a‏ 
9333330909 6 © 2 5 د دن 6 02 © © 6 © 5 56 6 / ER‏ 


که 


= | 


Ss — سم‎ Distributor ۳ 2.1 
Fiuidized. | ። 
31,388 4439 4 ات‎ ዘበ።ጋሪጋ ያያ ሪር የሪ 
Bed : LL ( 
Vall 





Teo Lest س م‎ መ لح تست سس‎ Rignt 
Vieu ) 








Front 


ገ ٔ Heat Zxcnanger Tune Design. 
All Dimensions ın ca. 





mixture of steel grit abrasive having a specific gravity of 8:1. 
Appendix E lists the size distribution of the steel grit used 
during the sampling operation. The bed was operated uithout 
recycling the fines captured in the cyclones. The static bed 
height of the material was 22.54 cm (8.875 in) throughout the data 


collection period. 
Equipment Set-ug 


Fig. 11 shows the placement of the sampling device inside 
the fiuidized bed. The sample trap was positioned such that it 
was directly over a spacing between tubes (22.86 cm (93 in) from 
the center of the froni-most tube] and-41.91 cm (16.5 in) from the 
left wall. The height of the sample trap above the distributor 
was varied during the sampling process as discussed in the section 


on sampling procedure. 


Fig. 12 shows the کے‎ of the bubble probe and the 
anemometer probe with respect to the sample trap. The bubble 
probe was placed directly below the sample trap and 26.57 cm (10.5 
in) above the distributor. The probe extension was placed at an 
angle so as not to interfere with the sample trap operation. To 
protect the anemometer wire from particles impacting it, a special 
shield consisting of * 320 mesh screen and an aluminum frame was 
slaced around it. ihe anemometer probe uas attached to the bubble 


probe extension with the entrance to the anemometer probe 29.21 cm 


(11.5 in) above the distributor. This placed the entrance to the 
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anemometer probe 2.54 cm (1.8 in) away from the center of the 
sample trap and 0.54 cm (0.25 in) outside the area directly below 
the sample trap. As a result, the probe has a minimum effect on 
the air flow from the bed to the sample trap, but will only 
measure the gas velocity at the edge of the sample trap perimeter. 
The positioning of the two probes above the distributor remained 


constant throughout the sampling evolution. 


Fig. 13 is a block diagram showing the equipment used during 
the sampling operations and their interconnections . Table l is a 
listing of the equipment used. The oscilloscope time base was set 
for MANUAL TRIGGER, SINGLE SWECP mode and a sweep time of SU 
ms/div. The channels of the cial trace amplifier were set at 5 
volts per division for the bubble probe and 2 volts per division 
for the sample trap inputs. On the differential. amplifier, one 
channel was not used and the second channel was set at 1 volt per 


division for the anemometer probe input. 
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Time 
Base 
Dual ies 
Trace Amplifier 
Amplifier 
Light 
Source 
Solenoid Optical 
Pauer Detector Anemometer 
Supply Power Supply 
/ኣ 
Bubble Anemometer Probe 
Probe 
ETS 3 Block Diagram of the Equipment Used 


During the Sampling Operations. 
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Oscilloscope 


Anemometer 


Optical Signal Detector 


TABLE 1 


TEKTRONIX 5111 Storage Oscilloscope 
TEKTRONIX SAÍ8N Dual Trace Ampl 
TEKTRONIX SAZ1N Differential Ampl 
TEKTRONIX 5218808 Time Base 


Thermal Systems Inc. 

1951-2 Monitor and Power Supply 
1954-A Linearized Anemometer 
Anemometer wire u/ #329 screen guard 


Optical Signal Detector Power Supply 


Optical Source and Power Supply 


Oscilloscope Camera 


Listing of equipment used during particle sampling 


operations. 
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To determine the fluidization conditions within the bed, a 
set of manometers were used. These manometers measured pressures 
within the bed, at heights from 4.13 cm (1.62 in) to 37.15 cm 
(14.62 in) above the distributor, in 2.54 cm (1 in) increments. 
The pressure data corresponding to each trap position and bed 


velocity is listed in Appendix I. 


a === i= To determine the gas flow conditions within the bed, an 
orfice fiow meter with 1D - 1/2D taps uas located upstream of the 
distributor. The computer program listed in Appendix G was used 


to convert the pressure tap data to mean air velocities within the 


bed. 


Sampling Procedure 


Data was collected for four (4) mean bed velocities at six 
(6) different sample heights. The sample trap was placed at a 
Given sampling height (measured from the distributor to the bottom 
of the trap), and ten (10) samples were collected at each desired 


velocity. The trap position was then changed to a new height. 


During certain sampling conditions, those which involved low 
sampling heights with the higher air velocities, the paddle arms 
would occasionally impact the sides of the sample trap and not 


close the sample trap completely. It is assumed that this occured 


ag 





when a large bubble erupted directly under the paddle arm and 
deflected the paddle arm into the side of the trap. Whenever this 


occured, the trap was de-energized and the closure cycle repeated. 


Each time the height or velocity was changed, a complete 
sampling cycle was conducted and this sample discarded. This was to 
prevent any accumulation of particles (in the entrance to the 
vacuum or purge lines on the sample trap) from being added to the 


first sample at the neu height or velocity. 


For each set of data at a given height and velocity the 


following information was recorded: 


1) Fluidized ር above the distributor 
determined visually and by pressure measurements. 

2) Pressure upstream of orfice plate (P1) 

3) Pressure difference across orfice plate (delta P) 

4) Air temperature in bed 

S) Pressure distribution in bed 


8) Height of sample trap above distributor 
The following procedure uas used during sampling: 
For each sample to be taken within a data set. 


Bein] .221l caneı tions 


1) Vacuum aır supply OFF 


2) Sweep trigger on oscilloscope RESET 
59 





3) 


4) 


B. Sampling 


1) 


2) 


3) 


4) 


O-scope memory ON-CLEAR 
Solenoid power supply switch 51 OFF 
Solenoid trigger switch $2 ON 
Procedure 


Trigger oscilloscope swees and wait until sweep 


is at the center of the CRT. 


Close the Solenoid power supply switch Si. When 
Si is shut, the oscilloscope will show an 
additional trace. This third trace is used to 
determine the closure time relative to the 
presence of gas jets and bubble eruptions. An 
example of a typical oscilloscope trace is shown 


ip Fig. l4. 


Open the solenoid trigger switch SZ. This 
reduces the current to the solenoid to prevent 
Overheating. The maximum allowed time to let SZ 


remain closed is five (5) seconds. 


Turn on the air supply to the vacuum system and 
leave cn one (1) minute. The exhaust air from 
the venturi on the vacuum system must be directed 
into the purge line in an oscillatory manner. 


This will agitate the particles within the trap 
51 
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so as to move them into the vacuum line. 
5) Turn off the vacuum air supply. 
6) Turn off the solenoid power supply switch Sl. 


7) Remove sample from sample container and place in 


specimen bottle. 
8) Photograph trace on oscilloscope. 


9) Repeat from A. 


Sample Analysis 
Sample Weight Determination 


Each sample obtained was weighed, using a Torsion Balance Co. 
TORBAL scale, to an accuracy of 0.01 grams. The average value and 
standard deviation uas then determined for each set (specific 
height and velocity) of samples. The ueight of particles in a 
completely filled sample trap uas also determined for void 


determination. These results are listed in Appendix I. 
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Determination of Particle Size Distribution 


To determine tne particle size distribution which occurs at 
various heights for a specified gas velocity, three (3) samples 
from each set of data at a specified height and velocity, and fron 
the bulk bed material were analyzed using a Zeiss Videoplan 2 
Image analyzer. The software used was the “Image Analysis System 


MOP-Videoplan" distributed by KONTRON Electronics Group. 


For each sample selected, a microscope slide uas prepared 


— 
— 


using double sided adhesive tape on uhich a portion of the 
selected sample uas placed. To ensure a sharp contrast was 
achieved, each slide was backed with white paper. The sample was 
then placed under a Zeiss microscope to which Hr image analyzer 
was connected via a vidiocon tube. The magnification used was 
SAX which provided an average view of about 8 particles at a 
time. The analyzer was then used to analyze the partices present 
on each slide of interest. The mode in which the image analyzer 
was used was the equivalent diameter mode. This mode determined 
the cross sectional area of each particle sampled and calculated 
the diameter of a circle with the same area. When a data set was 
completed, the data was analyzed for particle size distribution. 
Tne resulting output (Appendix J) consists of a particle count vs 
particle size histogram, a cumulative frecuency plot and a 


classification data list. A gaussian distribution fit for the 


data is also plotted on the histogram and cumulative frequency 


S4 





plots. The X axis of the plots are normalized with a range of 
zero (0) to four (4). To determine the actual diameter of the 
particle in microns for a given normalized value, the normalized 


value must be multiplied by the conversion factor 248.887. 
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COMPUTER MODEL 
Introd:c3ion 


This chapter descibes the theory, logic flow and testing of a 
particle trajectory computer model to predict the particle loading 
in the freeboard due to erupting buubles. The program, listed in 
Appendix L, is written in HP BASIC 2.0 and uas run on a Hewlett 
sa 9815 microcomputer. The output from the model is discussed 


in chapter VI and compared with experimental results. An in depth 


analysis of the program logic and structure is given in Appendix L. 


Model Theory 


The model developed here, is based on calculating the 
trajectory of a single particle as it is ejected from the bed 
surface and is acted upon by gravitational and drag forces. The 
drag force is due to the difference in absolute particle and air 
velocities. The air velocity is a combination of the initial jet 
velocity produced when a bubble bursts, and the superficial bed 
velocıtv. To ensure that ihe particle drag is calculated 
accurately, ihe follouing drag coefficient correlation given by 


White [17] was used. 
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Re 1+ YRe' (3) 


where: 
Cd = Orag coefficient for sphere 
Re Reynolds number 


Eqn 3 is valid over the range (Re 10. To calculate the 
particles position and velocity, the computer uses a forward 
difference method. Using Neuton’s Law (ZLF=ma) the acceleration of 


the particle due to gravity and drag is determined. Inserting 


-= 
- 


this acceleration into Eqn 4, the particles new velocity is 


determined. 


¥V=Vo tat (4) 
where: 
V = Particles new velocity 
Vo = Particles present velocity 
a = Acceleration of particle 
t = Time increment of calculation 


To determine the particles new position, the velocity 


eakulatced in Ecn 4 is inserted into Eqn 5. 
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H = Ho + (Vo + V) + 
ረ (5) 


where: 
H = Particles new height 
Ho = Particles present height 
Vo = Particles present velocity 
¥ = Particles newu velocity 
t = Time increment of calculation 


Using Eqns 3,4,and 5, the trajectory of the particle is calculated 
from the time it initially leaves the bed until the time that it 


returns to the bed. 


These calculations are repeated over a range of particle 
diameters from 80 to 570 microns. By determining the residence 
time of each particle within a specified height increment (ô H) 
above the surface of the bed, a particle 0.1 distribution 
above the bed is determined. The height increment (AH) used in 
the program is 2 cm. At the end of each time step when the height 
calculation (Eqn S) is completed, the counter representing the 
particular 2 cm height increment which the particle is in, is 
incremented by one. Each particle size has its own set of 


counters to allow indivicual particle analysis. 


The calculated density distribution is then weighted with the 


particle size distribution of the bulk bed material since the 
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probability of a given size particle being present at a specified 
heignt is dependent upon the number of particles within the 
system. This is accomplished by multiplying each height counter 
of a given particle size with the number of particles for that 
given size present in the input bed distribution. By summing the 
density values for each set of particle diameters at a given 
height over the entire freeboard of the bed, the overall particle 


density above the bed surface is determined. 


The model assumes that all the particles are ejected 
perpendicular to the surface of the bed and are initially ata 
uniform velocity. Because the model uses single particles for the 
analysis, the effects of multiparticle interactions are not 


included in the model. 


Testing of Program 


To evaluate the validity of the en tuo tests were 
conducted. The first test compared the height solution produced 
by the computer with a closed form solution. The second test 
involved running the program with different particle diameter and 
particle distribution height intervals to ensure that a valid 


sample size was being used. 
losed Form Solution 


To determine a closed form solution for particle height as a 
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function of initial particle velocity and superficial velocity Uo, 
a force balance was used. The forces acting on a particle are 
gravitational and drag. The gravitational force, Fl, is simply 
the volume of the particle multiplied by the particles density and 


the gravitational acceleration, and can be uritten as: 


5 (6) 


where: 


ዮ1 Force due to gravity 
Pp = Density of particle 
=. D = Diameter of particle 
= gravitational acceleration 


0 
| 


In order to get a closed form solution that did not involve 
non-linear differential equations, Stokes flow was used for the 
closed form solution only. The computer model used the Stokes 
equation only to compare results with the closed form solution, 
afterward, Eqn 3 was used. Using the Stokes drag coefficient 


relation, the drag force on a particle can be determined as: 


E 3 هلا‎ 2 Up) ሥ 0 (7) 


FZ = Force due to drag 

Uo = Superficial bed velocity 
Up = Velocity of particle 
Absolute viscosity of air 
= Diameter of particle 


oF 
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By inserting Eqns 8 and 7 into Neuton's Lau (I F=ma) and 


simplifying, the following differential equation is obtained: 


KEES CE (8) 


where: 


DRS 
ዘ 


Acceleration of particle 
= Velocity of particle 
- Ci = 182. 
Pp D^ 
C2 = C1 Uo - g 
boundary conditions: 


1) t=@ ቨ=ሀፀ 
2) t-0  X-9 


This second order linear differential equation can be 
solved using the given boundary conditions uith the resulting 


closed form solution given as: 
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CST O CITO =D 
El (9) 


where: 


2 
i 


Height of particle 


Using the initial conditions listed in table 2, the solution 
obtained using the computer model (maximum height= 15.270 cm, time 
to maximum height= 0.275 sec) was identical to three decimal 


places with the solution obtained using Eqn 3. 


' TABLE 2 
Superficial Velocity (Uo): 69.96 cm/s (2 ft/s) 
Initial particle velocity (Upo): 304.8 cm/s (10 ft/s) 
Particle diameter: 200 microns 
Particle density: 5000 kg/m? 
Time increment (computer): 0.001 sec 
Jet velocity: 0.0 cm/s 


List of parameters used to check computer 
calculations egainst closed form solution. 
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Sample Size Sensitivity Test 


To ensure that appropriate sample sizes were used to minimize 
errors due to coarse sampling intervals, tuo sensitivity tests 
uere run. One test involved changing the particle diameter 
interval from 18 microns to S microns. The second test changed 
the height sampling interval (AH) from 2 cm to Sem and then 1 cm. 
For each of the tests, the same initial conditions were input into 
the program. Table 3 shows the resulting output from the program 
listing the initial conditions and the resulting diameter versus 
پوس‎ height data calculated. Fig. 15 shous a plot of the 
calculated maximum height vs particle diameter data listed in 


table 2. Fig. 16 shows the bed particle distribution used for 


each of the tests. 


Fig. 17a shows the entrainment calculation using a diameter 
interval of 10 microns and a AH of 2 cm. Fig. 17b shows the 
same calculation using a diameter interval of 10 microns and a AH 
of 5 cm. The curve is not as smooth but still retains the same 
general shape. The peak of the curve shown in Fig. 17b occurs at 
a height of about 2@ cm whereas the peak in Fig. 17a occurs at 
about 18 cm. A semi-log plot of these curves would show that the 
slope of the line to the right of the peak would be larger for the 
data represented by Fig. 17b. The effect of maintaining OH at S 
cm but decreasing the diameter interval to 5 microns is shown in 


Fig. 17c. There 1s no readily detectable difference between Fig. 
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Listing of Input and Resulting Maximum 
Particle Heights uith Time to Maximum 
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the Increment Sensitivity Tests. 
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CEMO TIO. 175. In Fig. 17d, the effect of changing AH to 1 
cm and diameter spacing to S microns is shown. The difference 
between Fig. 17a and Fig. 17d is barely noticable and no 
detectable change in the slope to the right of the peaks is 
present. As a result of this analysis, the program uas operated 
with a AH of 2 cm and a diameter interval of 0 mierons. This 
reduced the calculation time to half of that required uhen using a 
diameter spacing or 5 microns when the same total diameter span 
was used and, as was seen in Figs. i7a and 17d, the difference in 


output dees not require the finer increment. 
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CHAPTER N 


EXPERIMENTAL RESULTS AND DISCUSSION 





Minimum Fluidization Velocity 


Fig. 18 is a plot of mean bed velocity (Uo) vs pressure drop 
through the bed (Pb). The velocity values were determined using 
the program in Appendix G and the pressure data listed in Appendix 
I. From this plot, the minimum fluidization velocity (Umf) for 
the bed conditions used during this study is determined to be 9.15 


m/s (8.5 Ft/sec). 


Entrainment Analysis 


. Table 4 lists the averaged sample ueights and their standard 
deviations for the samples (fippendix I) collected by the sampling 
apparatus. Looking at the standard deviation of the sample 
groups, the standard deviation is fairly large compared to the 
average values. However, visual observations of the fluidized bed 
in operation would suggest that a larger standard deviation would 
be expected. The shor: sample cycle (2 ms) and the bubble burst 
activity in the bed are the main reasons for this conclusion. The 


average standard deviation is 29 % of the average sample weight 
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and the range extends from 68 % of the average for the sample set 
numbers 196 thru 115 to only 7 X of the average for the samples 
285 thru 295. This variation is believed to be due to the nature 
of the bubble activity within the bed. It would appear that when 
a bubble erupts at a time just prior to closure of the trap, a 
sample weight greater than the average would be obtained. 
However, when no bubble has been present, the material caught by 
the trap should comsist mainly of particles returning to the bed 
and_ therefore be below the average ueight caught. Table 4 also 
shows the density of each sample collected. These values were 
determined by dividing the averaged sample weights by the volume 


of the sample trap which is 43.02 cubic cm. 


— 
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TABLE 4 


Superficial Bed Velocity (Uo): 58.1 cm/s (1.905 ft/s) 
Non-dimensionalised velocity (Uo/Umf): 3.81 
Sample Height Average Standard Average 
Numbers Above Sed Weight Deviation Density 
(cm) (grams) (grams) (grams/cm ) 
296-7358 تسه حا‎ (.229B-—-:c 0.51 1.68 (-@1) 
2465-255 © 4.48 8.72 1.94 (-01) 
1365-125 12.4 2.41 9.76 5.56 (-0ረ) 
165-175 17.8 1 8 0.15 24581 (0%) 
76-85 2 2 1232 0.23 3.07 (-02) 
176-185 کت‎ 0.31 6.10 7221-7 
Superficial Bed Velocity (Uo): 1۳8 2 cq/Owel.S585 ft/s) 
Non-dimensionalised velocity (Uo/Umf): 207 
Sample Height Average Standard Average 
Numbers Above Bed Weight Deviation Density 
(cm) (grams) (grams) (grams/cm ) 
۶ 75 - 5 SI A235 15 1.01 )- 01( 
236-245 Beg 1.88 0-27 ۸۰ 16 ۲ - 21 
96-195 13.6 1.36 9. 32 3.16 (-@2) 
1556-165 | 2፡1. 9.56 0.12 153930202) 
56-5 24.1 0.52 0.36 2.14 (-02) 
188-95 31.8 018 9.94 4.18 (-03) 
Superficial Bed Velocity (Uo): - 39.5 cms (1297 PMs) 
Non-dimensionalised velocity (Vo/Umf): 2053 
Sample Height Average Standard Average 
Numbers Above Bed Weight Deviation Density 
(cm) (grams) (grams) (grams/cm ) 
2EB-275 7 1.31 037 د3‎ 05 (202) 
ZUS- 285 55 0.52 ۵2-7 1.44 (-22) 
126-135 14.3 9.28 0 8 5 51 ع)‎ ۵ 
136-145 17 0 8 9.04 8620.93) 
8695 24.8 0. 25 0.05 25095005) 
195-265 32 4 0.06 0.2 1.39 )-03( 
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TABLE 4 (cont) 


Superficial Bed Velocity (Uo): 35.4 cem/is (1.161 ٦ 
Non-dimensionalised velocity (Uo/Umf): Zo E 
Sample Height Average Standard Average 
Numbers Above Bed Weight Deviation Density 
(cm) (grams) (grams) (grams/cm ) 
255-255 6.4 0.83 0.13 1.45 (-82) 
286-225 10.2 0.32 0 4 7.44 (-03) 
1106-7 14.3.5 0.221 سه ۷ېي‎ 05 5.11 (-03) 
146-155 20.3 0.02 . 0.01 4.65 (-@4} 
66-7S 26.0 07 0 4 2 - و9‎ 0-05) 
286-7215 3340 0.01 = 9.004 2.32 (-04) 


List of experimental data showing sample averages, standard 
deviations, heights, and velocity conditions measured. 
Density values are calculated by dividing the average samnle 
weight by the sample trap volume. 


Fig. 19 shows the relationship between the average particle 
density caucht in the sample trap and the sample trap height above 


the bed surface as a function of Uo/Umf. 


For the relatively lou fluidization velocities used during 
the data measurements, (maximum Uo/Umf = 3.81) the freeboard 
height can Be assumed infinite. Using this assumption implies 
that all of the particles return to the bed and none are 
elutriated, ie comolete reflux. Under these conditions, the 
following equation has been suggested by Lewis et al (313 and Kunii 
and Levenspiel [1] to model the particle loading within the 
freeboard. 
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Plot of Particle Density vs Freeboard Height 
as a Function of Uo/Umf. This data was 
Obtained Using the Particle Sampler. 
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Pr = Po exp{-a 1} (19) 


where: 
Pr = Particle density at heicht | 
Po = Particle density at (+)bed surface 
a = Characteristic particle decay length 
1 = Height above bed surface 


Table S shows the values for the parameters Pc and a obtained 


by linear regression for the curves shown in Fig. 13. 
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TABLE 5 


Uo/ Unf Po a Correlation 
Coefficient 


2 0.238 9.1097 02379 
ST 9.۶ 0.1023 0.344 
eres 0.048 0.1 02537 
2000 0 


.031 n 399 4 857 


Pre = Po expí-a 1} 
Resuits of linear regression analysis for particle 


loading density (grams/cm ) vs height above the bed 
surface (cm). 


The parameter Po physically represents the particle loading 
density which would be obtained if the sampie were taken at the 
surface of the bed. This is not necessarily the case as is 
indicated by the computer model which is discussed in chapter VI, 
but is oniy a parameter descibing the particle loading 
distribution in the region of the data obtained. The dashed lines 
below freeboard heights of 4 cm indicate the region in question. 
The parameter a is a characteristic length of decay for the 


particle flux. 


A correlation between ¿he values for Po in table 5 and Uc/Um: 
is shown in Fig. 20. This plot shows that Po is closely related 


with Uo/Umf. Po varies with (Uo/Umf - 1) approximately to the 2.9 
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power, where as Walsh et al [10] determined the coefficient to be 
about 2.1. Wen and Chen [5] reported a correlation for particle flux 
which 1s proportional to bubble diameter and (Uo/Umf - 1) to the 


5/2 power. These correlations are listed in table 6. 


TABLE 6 


== » iPresenti work: 8 := 8.-13-(ህ6/ህጸየ - 1) -- ፤ዩ6/ ሸጠ sec 
Wen and Chen B = 1.34 (04) (Uo-Umf ) kg/ m sec 
Walsh et al B = 18 (Uo - Umf) kg/ m sec 


Comparison of least square fit relations for Po 
as functions of Uo and Umf. 


The differences in these correlations are due to the 
different bed configurations in which the data was taken and the 
measurement technique used. The present work utilized a bed with 
a relatively closely spaced tube configuration and steel grit 
(S.G. 8.1, median size 230 microns) for the bed mass. Both 
ascending and descending particle fluxes were captured in the 
sample. The work of Walsh et al [10], used a bed with tuo (2) 
widely spaced horizontal serpentine tubes and Ottawa sand (5.6. 
2.6, median size 755 microns). Also, only descending particle 
flux was used in determining their relations. The correlations of 
Wen and Chen [5] are a result of studies conducted on previous 


research using cylindrical column beds and low mass bed materials 
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(S.G. 0.8 - 2.5). The data for these analysis is based mainly on 


pressure measurements. 


It has been observed by Lewis et al [3] and Wen and Chen [51] 
that a is not a strong function of Uo. Both of these studies 
recommend that the characteristic particle decay length, 1/a, 


could be approximated by an expression of the form: 


1/8 = C Uo (11) 


Table 7 is a list of correlations obtained by other studies 
and in the present work. The study by Lewis et al [3] was 


conducted with 7S micron glass spheres in a cylindrical bed. Fig. 


21 shows the relationship between 1/a and Uo in the present work. 


TABLES? 


Present work 1/2 < 40.19 + ۵.03 (0 m 
Lewis et al t/a = (1.42 + 8.14 s)Uo m 


Wen and Chen 1/a = (8.25 + 8.29 s)Uo m 


.05 s)Uo m 


+ 
ጩ 


Walsh et al l/a ۱0 ተ 


Comparison of less: scuare fit relations for 
1/a as functions ot Uo (m/s). 
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Particle Size Distribution 





Two methods for determining particle size distribution were 
used. The sieve method, described in Appendix E, was used for bed 
material analysis only, due to the small sample sizes obtaincd from 
the trac. To analyze the small trap samples, an image analyzer 
was utilized. The procedure used with the image analyzer is 


described in the experimental procedure chapter. 


The sieve data representing bed mass as a function of 
particle size, is listed in Appendix E. Fig. ZZa shows the bed 
mass data converted to particle number as a function of particle 
diameter where the interval of particle diameter between 
successive measurements is 10 microns. These plots are faired 
from the data listed in Appendicies E and J. Fig. 22b shows the 
bed mass distribution. The conversion from a mass distribution to 
a particle number distribution was calculated assuming that each 
particle was spherical in shape. The volume corresponding to a 
given particle diameter was multiplied by the particle density to 
get a unit particle mass. The mass fraction of the sieve analysis 
corresponding to the specified particle diameter was then divided 
by the unit particle mass to obtain the representitive particle 
number. The overall resulting particle number distribution curve 


uas then normalized uith rescect to a maximum value or 190. 


The particle size distribution data obtained from the image 
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analyzer is listed in Appendix J. The image analyzer was used to 
analyze the bed mass and sample data from the Uo/Umf = 3.81 data 
set. The bar graphs and the cumulative percentage plots show the 
number of particles viewed by the image analyzer plotted as a 


function of normalized particle diameter. 


Fig. 22a shows a plot of the particle number distribution as 
a function of particle diameter as determined by the image 
analyzer. The data is normalized with respect to 100 and is 
compared with the data as determined by the sieve analysis. The 
descrepancy between the two piots can be explained by the methcds 
used to determine the respective data. The analyzer firs? 
cetermines the cross sectional area of the viewed particle. A 
circle, having -the same cross sectional area, is then calculated. 
_ This results in averaging the smaller minimum diameter with the 
larger maximum diameter of all particles with long cylindrical or 
ellipsoidal shapes. In the sieve however, a large portion of the 
particles will pass through the sieve screen By means of the small 
cross sectional area presented by their longitudinai direction. 
To be consistant in the following sections, the image analyzer 


data wiil be used to correlate all particle density distributions. 


The particle size distribution as a function of height above 
the bed surface was evaluated using the data collected at Uo/Umf = 
3.81. The image analyser ues used to determine the distribution 


by vieuing three (3) random samples from each of the six (B) 


height positions. Table 8 lists the various statistical values 
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tor the particle size distribution data obtained (Appendix J). 


TABLE 8 
Freeboard Average Median Mode Diameter at 
Height Diameter Diameter Diameter SOX of 
Cumulative 
Number 
Distribution 
( cm) (microns) (microns) (microns) (microns) 
3,8 193 192 142 153 
7.6 219 214 187 203 
123 204 189 162 157 
18 173 164 154 179 
22.7 187 179 152 I 
|1 212 137 134 139 


Staistical values for particle number distribution 
as a function of freeboard height. A complete 
listing of the data is given in Appendix J. 


Fig. 23 shows these values plotted against the bed height at 
which they were taken. The excursion of points at the 7.8 cm 
height is assumed to be due to analysis error. The linear 
regression lines for the average and median values show that they 
are weak functions of collection height. For the particle 
diameter values representing the mode and the S@Z point on the 
cumulative percentage plot, the linear regression lines show a 


stronger dependence on collec:zion height. 
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Oscilloscope Trace Analysis 


Figs. 24 thru 27 show examples of the oscilloscoge traces 
given in Appendix K. Figs. 24 and 25 represent the traces 
obtained during low Uo coditions. Under these conditions, the 
normal bubble probe output is high (+ 5 volts dc). When a bubble 
errupts, particles are throun from the surface of the bed and 
eclipse the light path at the tip of the probe (trace goes to 
zero). Figs. 26 and 27 represent the traces obtained at higher 
Uo conditions when the bubble probe is normally eclissed (bubble 
probe output is lou, zero) by particles. The presence of a bubble 
is detected by a high output from the bubble probe due to the 
bubble creating a void through which the light beam can pass. The 
anemometer output just above the bubble probe is used to determine 
the velocity of the gas jet leaving the erupting bubble. As seen 
in the figures, there is not always a bubble associated with a gas 
jet and vice versa. This is due to the anemometer probe being 


Slightly offset from the bubble probe as described in the 


experimental procedure section. 


The gas jet in Fig. 24 is delayed S ms and the gas jet in 
Fig. 25 is delayed 35 ms from the point where the bubble traces 
begin their excursion to the zero (eclipsed) condition. Since the 
seperation between the bubble probe and the anemometer is 1 œm, 


this results in an estimated jet velocity of 200 cm/s (Fig. 24) 
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and 28 cm/s (Fig. 25). Using the calibration curve in Appendix 
M, the velocities are determined to be S18 cm/s for the jet in 
Fig. 24 and 792 cm/s for the jet in Fig. 25. The large 
difference in these values could be due to a bulge in the bed 
surface from the bubble underneath ecligsing the bubble probe. 

If this is the case, there would be a delay between the time at 
which the bubble probe wes eclipsed and the eruption of the bubble 


from the bed surface. 


An analysis of Figs. 26 and 27 results in the same 
discrepancies betueen velocities calculated from delay times and 
measured by the anemometer crobe. Fig. 27 shows another phenomena 
which occurs quite frequently. The gas jet appears before the 
- bubble is detected. : There ود‎ several possibilities which can 
explain these occurances. First, the bubble occurs off center 
from the bubble probe. Under this condition, the bubble may erupt 
and initiate a et which is registered by the anemometer. The 
bubble then continues to rise and the side of the bubble is 
registered by the bubble probe. This explanation can be altered 
to include bubbles coalescing belou the surface. Under this 


condition, a bubble may be pulled into the vortex of an already 


eruoting bubble and it eclipses the bubble probe. 


An analysis of the jet velocities:and durations for all the 
oscilloscooe traces in fippendix K, resulted in the cata listec in 
tai 3 j t lues to 5 d for the jet 
table 2. In later calculations, the values to De use t j 


velocity and duration will be 609.6 cm/s (20 ft/s) and 20 ms 
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respectively. 





TABLE 9 
Average Stnd Dev 
Jet Velocity 689 cm/s 559 cm/s 
22.5 ft/s 18. ft/s 
Jet curation 21.25 ms 13.22 ms 


Average and standard deviation of jet velocity 
determined from oscilloscope traces in 
Appendix K. 


Versus Bed Sctivityv Correlation 


The main purpose of inserting the bubble and anemometer probe 


- Beneath the 


correlation 
presence of 


photographs 


sampling apparatus- was to determine whether a 
exists between the sample weight collected and the 
bubble eruptions and gas jets. To evaluate the 


taken of the oscilloscope traces (Appendix K), the 


following information was required: 


A. The average velocity of the particles as they 


travel from the bed surface to the trap. 


B. The distance the particle must travel to reach 


the trap. 


C. The closure time of the trap relative to the 
bubble eruptions and gas jets. 
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The average particle velocity was calculated by determining 
the average particle size and using the height output from the 
computer trajectory model. The average narticle size as 
determined by the image analyzer uas 190 micons. To determine the 


average velocity using the model, the 


follouing initial conditions 


1 


vere input to the progran: 


superficial velocity (Uo) 57.9 m/s (1.4 ft/s) 


This was the actual velocity 
measured Guring samoling. 


Jet velocity (Uj) 639.5 و وی‎ (20 ft/s) 


This value was determined from 
the anemometer data. 


Jeb duration (t jet) 22 ms 


This value was determined from 
the oscilloscope traces. 


The height attained by a 180 micron particle, as determined 
by the model with the above conditions is 44.7 cm in @.295 
seconds. This results in an average velocity of 44.7/0.285 - 


151.5 om/s (4.97 ft/s). 


1۳۲2 Gistame= a perticl= wust travel to reach the center of 
ሽብ wees trom the bec surface is cbteined from the cata listed in 


Resendıx I. Since trag height is measured from the bed surface to 


the bottom of the trap, 4 cm (9.13 ft) must be added “o the trac 
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heights to obtain the distance to the center of the trap. The 
total time required for the particles to leave the bed surface and 
arrive at the center of the trap can now be determined. Table 10 
shows the transit time for each set of data for which the 


oscilloscope traces were photographed (Appendix K). 


TABLE 1@ 

Data Trap Trap Total Time 
set Height Height Transit Before 
Number Bottom Center Time Actuation 

(cn) (cm) (ms) (ms) 
56-65 24.13 22.13 173 136 
66-75 26.03 30.93 198 155 
25-85 2 23 2623 173 130 
86-35 24.76 28.76 190 147 
96-105 13865 122605 11) i 3 


List of transit times for particles traveling from 
the bed surface to the center of the trap. The time 


prior to actuation of the sample trap is also shoun 
(Total time - 42.6 ms). 


The total closure time of the sampling apparatus (closure of 
solenoid power supply switch SI until trap is shut) is determined 
in Appendix D to be 42.6 ms. This value is subtracted from the 
garticle transi* time to obtain the time before actuation in which 


particles leaving the bed surface will be caught in the trap. 


To determine whether or not a correlation exists between the 
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weight of a sample and the amount of bed activity present in the 
bed, a weighted analysis was used. Each sample was evaluated 
three ways. First, if the sample weight was less than the average 
sample weight, a weight factor (U) of -1 was assigned to it. If 
tne sample weight was greater than the average, a weight factor of 
*l uas assigned. Second, the oscilloscope trace was analyzed at 
ihe point corresponding to the "time before actuation” listed in 
Table 18. This position is the time before the solenoid power 
supply trace is present (Fig. 28). If the trace showed signs of 
a bubble eruption or a gas jet at this point, the position factor 
(p) uas assigned a value of +l. If no activity was present, the 
position factor was assigned a value of -1. Third, each trace was 
analyzed again at the “before actuation time" but bed activity 
within a +/- 2S ms region was counted. If there was bed activity 
within this region, the region factor (r) was assigned the value 
of +1. If no activity was present, the region factor was assigned 


the value of -1. The following equation was used to assign a 


correlation factor to each sample: 


Q=Up+Ur (10) 
where: 
Q = Correlation factor 
W = +1 if sample is ) average of data set 
-1 if sample is ( average of data set 
p = -1 if no activity is present in bed 


+1 if bubble or gas jet activity present 
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r = -i if no activity is present +/- 25 ms 
+l if activity is present +/- 25 ms 


Table 11 shows the correlation factor for each sample. The 
resulting average correlation factor for all samples is 0.32. 
This value suggests that no correlation can be made between the 
sample weight and bed activity. There are several explanations 
for this. Bubbles not detected by the bubble probe or debris from 
previous bubbles returning to the bed surface can influence the 
sample weight by increasing the amount cf particles collected. 
Also, bubbles of smaller size have a smaller velocity and would 
therefore produce a particle stream which is either slower or just 
delayed in leaving the bed surface. This would result in lower 


sample weights than expected. 
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TABLE 11 


Sample 

Number لما‎ a E Q 
55 1 1 1 2 
27 a 1 1 -2 
58 1 1 1 d 
c3 -1 -1 -1 2 
sa 1 -1 -1 -2 
=ረ Een 1 1. =>. 
53 1 1 1 ረ 
64 p -1 -1 Z 
56 E -1 1 0 
67 1 1 1 2 
69 E 1 1 -2 
78 = 1 1 2 
7 1 1 1 2 
ms >=! 1 ! -2 
83 1 1 1 Z 
84 1 1 1 ረ 
86 - 1 -1 1 0 
88 1 1 1 É 
83 1 1 1 -2 
20 21 ار‎ 1 0 
52 = -1 1 0 
55 1 1 1 2 
34 x. =1 1 0 
35 1 1 1 2 
96 t -1 1 3 
357 -1 -1 -1 ረ 
98 1 zi 1 Q 
Sg 1 =] =i ሙሪ 
109 1 =1 er m^ 
191 | 1 1 2 
192 SL = 1 1 3 
Average = ۵ئ‎ 

List of samples and their correlation parameters. The 


resulting average value for Q indicates that no 
correlation can be made from the data obtained to 
indicate by sample weight whether or not any bed 
activity occured below the sample trap. 
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CHAPTER VI 


TRAJECTORY MODEL RESULTS AND DISCUSSION 
Selection of Baseline Parameters 


The input to the model consists of the following five (5) 


parameters: 


1) Superficial bed velocity (Uo) 
~. 2) Initial particle E ( Upo? 

3) Peak gas jet velocity (Uj) 

4) Gas jet duration (tj) 


5) Particle distribution of the bed mass 


--- The Baseline values for each of these inputs uas determined 
to be as close to the actual experimental bed conditions as 


possible. 
Superficial Bed Velocity 


The superficial bed velocity was determined directly from the 
experimental data. For the bed conditions discussed in this 
section, the superficial velocity used is the same as the highest 
velocity condition under which the particle sampler was used. 


The superficial velocity within the bed for these samples was 
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calculated to be S7.9 cm/s (1.9 ft/s). 


Initial Particle Velocity 


To determine the initial velocity of the particle, it was 
assumed that the particle was located at the nose of a bubble and 
would therefore have the bubbles velocity. To determine the 
buböle velocity, Kunii and Levenspiel [1] give the following 


equation: 


y. 
Ub = 0.7M tg UBI EI عون د ون‎ (11) 


where: 
Ub = Bubble velocity 
Uo = Superficial bed velocity 
2 ro So ሀጠየ تت‎ Minimum fluidization velocity 


: . Ob = Bubble diameter 
SSS SS SES CS -ሀ graviatational acceleration 


Observations of the bed material during the sampling 
operation suggests that the average bubble diameter present in the 
bed is approximately 6 cm (2.4 in). For Uo = 57.9 cm/s, Umf = 
15.2 cm/s and Ob = 6cm, the bubble velocity is calculated to be 
97.2 cm/s (3.19 ft/s). This value was used as the initial 


particle velocity for the base line data. 
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Peak Gas Jet Velocity 


To determine an average peak jet velocity, the anemometer 
output on the oscilloscope traces (Appendix K) were analyzed. To 
determine the relationship between the voltage output from the 
anemometer (which is displayed on the oscilloscope) and the actuai 
gas jet velocity, the anemometer was calibrated in a wind tunnel. 
The resulting calibration curve for the anemometer output, as 
displayed on the oscilloscope trace, is given in Appendix M. The 
average peak velocity of the anemometer traces analyzed in chapter 


V was determined to be approximately 608.5 cm/s (20 ft/s). 


— 
— 


Gas Jet Duration 


S. To-determine an average gas-jet duration time to use as a 


baseline, the anemometer traces in Appendix K were used. The time 
duration of each gas jet analyzed is determined directly from the 
oscilloscope trace. The average duration of a gas jet uas 


determined to be approximately 20 ms in chapter V. 


Particle Distribution of the Bed Mass 


Two particle distributions for the bed mass were available 
for use as the baseline data. Both sarticle distributions were 
determined using the same sample but analyzed using different 
analysis techniques. The two methods used are the sieve method 


and the image analyzer method, both of which are described earlier 
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in chapter V. The two distributions are shoun in ۲ 7 
The bed distribution determined by the image analyzer uas used to 


be consistent uith the particle distributions measured at various 


freeboard heights. 


Tvpical Output Using Baseline Parameters 


Fig. 30 shows the resulting maximum particle height obtained 
as a function of particle diameter for the baseline conditions. 
The general shape of the curve is determined by the dominating 
force acting on the particle. For 2 particles, the dominating 
force at these air velocities is gravity. Therefore, momentum 
(initial particle velocity) is the controlling factor determining 
the maximum height attained by the particle. As the particle size 
decreases, the proportion of drag force to gravitational force 
becomes larger. For the smaller particles, the drag force becomes 
the dominate force. The result is that the maximum height 
attained by a particle contirwes to increase as particle size 
decreases. This continues until the value of Uo approaches the 
terminal velocity of the smaller particles. At this point, the 
particle is totally dominated by the air flow in the bed. This is 
seen in the small decrease in height attained by particles less 
than 118 microns as a result of the particles rapid deacceleration 
tc Uo after the gas jet has stopped while the larger particles 
continue a little higher cue to momentum. The particle height 


continues to decrease for a short time until Uo becomes equal to 
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Faired Data From Appendix E and J. Interval of 10 microns. 
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or greater than the terminal velocity for the remaining particles. 


These particles would be elutriated from the system. 


Figs. 31 thru 36 show the individual particle densities at 
difrerent heights above the bed. All of the heights shou an 
increase in particle density for the smaller diameter particles. 
This is the result of the smaller particles falling at their 
respective terminal velocities which is slower than for the larger 
particles. ihe result is that the smaller the particle velocity, 
the longer the particle exists uithin a given height region and 
+ 


thus, the larger the particle density. It is also due to the 


larger number of smaller particles present within the system. 


Fig. 37 shows the particle density distribution above the 
bed. The two curves represent the effect of varying the range of 
particle diameters used in the model. The two ranges are 8@ - S7@ 
microns and 80 - 1070 microns. The peak of the 80 - 570 micron 
distribution is at about 18 cm of bed height uhereas the other 
distribution peaks at about Lil cm. The region shown on the graph 
belou these heights is the area knoun as the splash zone. The 
shape of this part of the curve is due to the initial acceleration 
of the particles leaving the surface of the bed by the gas jet. 
The particles then begin to slow down due to drag and gravity 
resulting in an increase in particle density. The peak on this 
curve coincides with the maximum height attaıned by the largest 
particles, and therefore, it can be assumed that the change in 


slope is due to the loss of particles as they return to the bed. 
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Fig. 38 is a semi-ln plot of the same data as plotted in Fig. 

37. The sharp drop at the right end of the plot is due to 
analyzing particles with a minimum diameter of 80 microns. Had 
smaller particles been included, the slope would have approached 
zero instead of infinity which would model the elutriation of 
particles. The slope of the line to the right of the peaks shown 
in Fig. 38, which are located at 1! and 18 cm of height above the 
bed surface are the same and are approximately -0.117 grams/cm . 
To decrease computation time, the 80 - 570 micron particle 
distribution uas used in the follouing sensitivity analysis as the 


decay slopes are equal. 


Model Sensitivity Analysis 


222 2337 An analysis of the change in bed characteristics due to a 


change in one of the five-model input’ parameters was conducted. 


The parameters changed and the values used are as follows: 


Superficial velocity (Uo): 


Baseline Uo: 57.9 cm/s (1.9 ft/s) 


1. 39.5 cm/s (1 ft/s) 


2. 91.4 cm/s (3 ft/s) 


Initial particle velocity (Upo): 


Baseline Upo: 1 2 2 ems (3.13 ft/s) 


t. 61.0 cm/s (2 ft/s) 


AS Sos 
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3. 305.0 cm/s (10 ft/s) 


Gas Jet velocity (Uj): 


Baseline Uj: 609.6 cm/s (29 ft/s) 


1. 305.0 cm/s (10 ft/s) 


2. 457.2 cm/s (15 ft/s) 


Gas jet Guration (tj): 


Baseline tj: 20 ms 
1. 10 ms 
2. 30 ms 


Particle distribution of bed mass: 


. Baseline distribution: Image Analyzer 


7 E 1. Sieve analysis data 


The bed characteristics which are evaluated: against the 


baseline characteristics are: 


The maximum height attained by each particle from 80 to 


S70 microns. 


The density of each particle size (80 - 570 microns) at 


freeboard heights of: 


a. 4 cm (1.S in) 


0 ت 8 
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c. 12 cm (4.7 in) 
d. 18 cm (7.l in) 
e. ZZ cm (8.7 in) 


. Bi em (IZ Zin) 


3. The particle density distribution in the freeboard. 


Variation of Superficial Velocity (Uo) 


Fig. 39 shows the influence of Uo on the maximum height 
attained by the particles. As was described in the previous 
9ھ‎ the larger particles are dominated by momentum and not 
drag. _This is shown in Fig. 32 by the small change in maximum 
- height attained by the large particles due to a change in Uo. As 
. the particle size decreases, the drag force becomes the dominate 
force and the effect of Uo on particle height increases. As seen 


in Fig. 39, a change in Uo produces a moderate change in the 


maximum height attained by the small particles. 


Figs. 40 thru 45 shou the effect of Uo on individual 
particle densities at different heights above the bed. All of the 
heights show an increase in particle density for the smaller 


diameter particles as Uo is increased. This is the result of the 
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smaller particles falling at their respective relative terminal 
velocities which is slower for smaller particles. A decrease in 


the maximum particle diameter present at a given freeboard height 


as Uo is decreased can be observed in Figs. 4@ thru 4S. This 


198 





Freeboard Height (cm) 





| د7‎ cm/s 
| 


















400 500 600 


Particle Dia. (Microns) 


ہے << 


100 zoo 366 


| 
ro — و‎ 
| 








— س‎ a- 


8569 Maximum Particle Height vs Particle 


NE omm ٦ Diameter as a Function of Uo. 


Relative Particle Number 


> eS 
ہہ وچ‎ 27.9 cm/s 


| 
| 
| A 
| 
| 








ሻ 1 ቋር? ل 300 لالات‎ Sud 3 
Particle Dia. (Microns) 
Fig. 29 Relative Particle Number vs Particle 


Diameter as a Function of Uo. 
Freeboard Height of 4 cm. 


199 





ር. 
ር) 


Rel 


Relative Particle Number 


aj So om Sreeboard Height of 8 cm. 


ative Particle Number 


4 
لٹ‎ © 
E ——_ AMAN 4  :- SY AAA 


91.4 cm/s 


97.9 cm/s 
77771 30.5 cm/s 







3 100 "260 360 3180 ۰ 55 ۰ 


Particle Dia. (Microns) 





Fig. 41 Relative Particle Number vs Particle 


Diameter as a Function of Uo. 





120 

38 

22 91.4 cm/s 
[^ We 96 
zo بس چک‎ "771 30.5 5 
5 













12:2 2058 28 د‎ S^ 0 
> 


Particle Dia. (Microns) 
ES 2 Relative Partıcle Number vs Particle 


Diameter as a Function of Uo. 
Freeboard Height of 12 cm. 





Relative Particle Number 


Relative Particle Numbe 


)0 
لكا 


و یہ 
ወጪ‏ © 
A‏ —— 


1 سا‎ ከነ) ርህ Ab N 
ፕ መመ መ O O لک‎ 


a. 


En w3 5 


108 





.  —Freeboard Height of 18 cm. 


91.4 cm/s 


=>" 57.9 cms 
— 30.5 cm/s 








IHD 200 SG 400 Se 052/27 
Particle Dia. (Microns) 


Relative Particle Number vs Particle 
Diameter as a Function of Uo. 





مج مہم nn‏ — 





cm/s‏ 91.4 سب 
cm/s‏ 57.9 ہم ہے — 





1409 ¿BO 28393 400 0 BUN 
Particle Dia. (Microns) 
Relative Particle Number vs Particle 


Diameter as a Function of Uo. 
Freeboard Height of 22 cm. 





91.4 cm/s 


e O TE 
-- - - - س‎ 38.5 cm/s 


Relative Particle Number 





Y ۱9 200 300 300 Sug 9 
Particle Dia. (Microns) 
Fig. 45 Relative Particle Number vs Particle 


Diameter as a Function of Uo. 
Freeboard Height of 31 cm. 





decrease is due to the decrease in lift given to the ascending 
particle as Uo is decreased. As a result, all particles achieve a 


lower maximum height when Uo is decreased. 


Fig. 46 shows the effect of Uo on the particle density 
distribution above the bed. There is a slight increase in overall 
narticle density with an increase in Uo. As expected, there is 
also an increase in the maximum height attained by the particles 
with an increase in Uo. Fig. 47 is a semi-ln plot of the same 
data plotted in Fig. 46. Table 12 lists the slopes of the lines 
to the right of the peak shown in Fig. 47 which is located at 
about 18 cm of height above the bed surface. This table shows an 


^ 


increase in the particle distribution slope as Uo is decreased. 


TABWE 12 
Uo Slope Intercept 
cm/s gms/cm gms/cm 
30.5 -0.136 7.036 
5 د2‎ zu 4 6.758 
91.4 -0.098 6.473 


Effect of Uo on the slope of the 
particle density distribution as a 
function of height for the 
distributions shown in Fig. 47. 
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Variation of Initial Particle Velocity (Upo) 


Fig. 48 shows the relation between initial particie velocity 
(Upo) and the maximum height attained by particles of different 
size. This figure also shows that the larger particles are mainly 
momentum dependent while the smaller particles are drag dependent. 
If Upo were increased further, the larger particles would continue 
to increase their maximum height. The smaller particles would 
approach a height limit wnich is dependent upon the superficial 


velocity in the bed. 


Figs. 49 thru S4 show the effects of increasing Upo on the 
individual particle densities for increasing heights. In general, 
the individual particle distributions undergo the same relative 
changes from the bed surface to the maximum height position. The 
difference m the height above the bed surface at uhich the 


particular distribution is present. 


Fig. 55 shous the effect of varying Ugo on the particle 
density distribution in the freeboard. fis uas noted previously, 
the peak density occurs at the point uhere the largest particles 
attain their maximum height above the bed. As Upo is increased, 
the curve to the right of the peak becomes shorter and steeper. 
This trend will continue until the the effect of the smaller 


particles returning to the bed cause the slope to decrease. Fig: 
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58 is a semi-1n plot of the data in Fig. 55. Table 13 confirms 
that the slope is becoming steeper as Upo is increased although 
the influence of initial particle velocity on the slope of the 


density vs height curve is moderate. 


TABLE 13 
Upo Slope Intercept 
cm/s gms/cn gms/cm 
61.9 =A 118 6.958 
gura -0.117 6.758 
152.4 -0.125 8.500 
304.8 -0.299 18.20 


Effect of Upo on the slope of the 
— سس‎ S particle density distribution as a 
لل‎ m “function of height for the 

distributions shown in የ16. 58. 


Variation of Jet Velocity (Uj) 


Fig. 57 shows the effect on maximum particle height when Uj 
is varied. As in the two previous analysis’, the influence of 
momentum and drag on the different particle sizes is apparent. A 
decrease in jet velocity from the baseline value results ina 
large drop in small particle height and a small change in large 
partıcle heıcht. This suggests that jet velocity has a strong 
effect on the Giscersion, or seperation of particles of different 
sizes at increasing freeboard heights. Increasing Uj results ina 


stretching effect of the particle size distributions above the 
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bed. One interesting point to note is the loss of the rounded 
edge at the heights corresponding to the smaller particles. This 
shous the extremely large effect that drag plays when the air 


velocity is very close to the terminal velocity of the particle. 


Figs. 58 thru &3 show the effect of varying Uj on the 
individual partical densities at varying heights above the bed. 
Again, as in the previous parameter analysis, the same 
distribution shapes can be seen for each value of Uj, the only 


difference being the height at which it occurs. 


~ Fig. 64 shows the change in particle density distribution as 
= == -Uj.is-varied.--rig. =65-is-a-semi-In-plot of this data showing 


that as Uj is decreased, the slope of the distribution changes 


substantially. Table 14 also shows this effect.: 


TABLE 14 
Uj Slope Intercept 
cm/s gms/cm gms/cm 
503.5 -0.117 6.758 
45702 -0.207 7.380 
304.8 - 0.400 7.300 


Effect of Uj on the slope of the 
serticle density distribution as a 
function of height for the 
distributions shown in Fig. 55. 
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Variation of Jet Duration (tj) 


Fig. 66 is a graph showing the effect of varying jet 
Guration on the maximum heights attained by the particles. It is 
easily seen that a small change in jet duration results in a 
substantial change in particle height. This effect indicates that 
the distribution present in an actual fluidized bed is probably a 
statistical average of a rapidly fluctuating particle distribution 
wnich is controled by the durations of the jets from the 


neighboring bubbie eruptions as well as local eruptions. 


- 8 — መመጮ፦ 


Figs. _67 thru 72 show the effects on the individual particle 
density distributions, at different freeboard heights, as a 
function of jet duration. These figures also show the drastic 


change in density distributions caused by changes in jet duration. 


Fig. 73 shows the particle density distributions above the 
bed as affected by jet duration. Fig. 74, which is a semi-In 
plot of the data, shows the drastic changes in the decay slopes. 


Table 15 also shows this drastic variation in slope. 
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TABLE 15 


ti Slope Intercept 
cm/s 3.5 / ርሾ 185 / 5 
19.9 پت‎ 12.0 
20.9 0 6.758 
29.6 -3.020 6.23 


Effect of tj on the slope of the 
particle density distribution as a 
function of height for the 
distributions shoun in Fig. 74. 


luu Variation or Particle-Distribution in Bed Mass 


vues 31:6 cue rige (Seshous=the tuo-bed. particle distributions used in 


ፒዑ. ጨ--“ o — = 


this analysis. The determination of these distributions is 


8 $e descibed earlier in this uork. 


Figs. 


76 thru 81 show the effects of varying the bed 


particle distribution on the individual particle distributions at 


different heights above the bed surface. Using the sieve particle 


distribution instead of the image analysis distribution produces a 


shift in the mean particle distribution towards the smaller 


perticles. 


The distributions resulting from the sieve data do not 


shift as much as the image analyzer data uhen the freeboard height 


is increased. The small shift in mean diameter exhibited by these 


pauses coWsistan: Wicm the distributions listed in Appendix J. 


F10- 


82 shows the particle density distributions in the 
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freeboard for the tuo bed mass conditions. The effect of the 
Steve distribution (which has a much lower average particle 
diameter) is to decrease the initial rate at which particle mass 
returns to the bed. This is due to the smaller particles having a 
mich smaller terminal velocity. This will not only raise the bulk 
of the particles to a higher height, but will also increase the 


time required for the particles to return to the bed. 


Comparison of Model with Experimenta! Results 


As discussed in the previous section, the baseline parameters 


— 


— 


were selected as the experimental conditions present in the bed 
when sampling at the highest Uo setting. Table 16 lists these 


- Baseline - parameters again for review. 


TABLE 16 


Uo 57.9 cm/s (1.9 ft/s) 
Upo 97.2 cm/s (3:19 ft/s) 
Uj = 609.6 cm/s (20 ft/s) 
tj 20 ms 


Baseline parameters used in computer model. 


Figs $83 shows he particle density distribution predicted by 
the model using the baseline conditions. Fig. 384 is a semi-log 


plot of this same data. A comparison of the slope of the 
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distribution line after the peak value in Fig. 84 with the slope 
obtained from the experimental data shows very good agreement 
between the two. This data is also given in Table 17 along with 
the cata obtained for the other experimental data sets and is 
shown plotted in Fig. 85. The slopes for the model data show a 
steady increase in slope as Uo/Umf is decreased. This trend is 
also followed by the experimental data but is not as smooth. The 
slope of the experimental data is aiso changing faster than the 
slope of the model data. One reason for this is that all the model 
calculations were performed using the same value for jet velocity 
and duration of 609.6 cm/s and 2@ ms respectively. The magnitudes 
ot_these values should decrease with decreaseing Uo/Umt, t the 
values to be used for these different concitions could not be 


determined. 


TABLE 17 

Uo/Umf Experimental Model Diff 
Slope Slope y 

gms/cm gms/cm 
دو‎ -0.1097 0-117 6.6 
dM T -0.1029 -9.173 19.5 
2 9 zu. ٦1 -0.129 ۹.7 
202 -0.1399 -0. 1 32 55 


Comzarison of slopes for the particle density 
distributions 2cove the bed as derived from 
the experimeniai cata and the computer model. 
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EG oo Comparison of slopes for the particle 


density distribution above the bed 
as derived for experimental data and 
computer model output. 
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The slope shown in Fig. 84 for heights between approximately 
4 and 18 cm is totally different from the slopes listed in table 
17. This is also the lower range in which the experimental data 
was collected. As was shown earlier, when the particle 
distribution analyzed included particles from 80 - 1070 microns, 
the peak was located at 11 cm. This is the height uhere the 
largest particles attain their maximum height and begin to fall 
back to the bed.. It does not seem likely that including the feu 
particles present above 1070 microns would totally account for 
this discrepancy. Other factors which the model has not accounted 
for wrich may explain these discrepancies is the effect of varying 
jet- velocities and durations.uhich were shown in chapter V to vary 
Greatly, and the effect of particles whose velocity vector is not 


Perpindicular to the bed surface. 
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CHAPTER VII 


CONCLUSTONS AND RECOMMENDATIONS 


Conclusions 


Several important conclusions can be made concerning the 


particle trajectory computer model and the experimental data. 


— 1. A particle trajectory model, even as simple as the one 
discussed in this uork, produces results uhich closely predict 
several aspects of particle activity within a fluidized bed. The 
more important predictions include the particle density 
distribution in the freeboard and the height distribution of 
particles. The results obtained with the current model indicate 
that further work on improving the particle trajectory model to 
include particle-particle interactions and paticle velocities not 


perpendicular to the bed surface is very desirable. 


2. The sensitivity analysis indicates that the jet duration 
and jet velocity are critical parameters in determining particle 
loading conditions in the freeboard. Since these are constantly 
changing from one bubble eruption to the next, a statistical 
distribution will be required to accuarately model freeboard 


particle activity. The developement of a basic model to predict 
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the jet Guration and velocity as a function of bubble size, Uo, Umf, 


etc is therefore needed. 


3. Particle distribution analysis of experimental samples at 
increasing heights above the bed, show the presence of large 
particies. Many of these are above the maximum calculated 
trajectory heights for these particle sizes. This indicates that 
either the jet velocity or jet duration carried these particles to 
these heights, or, particle-particle collisions are present in 


numbers great enough to be important in the analysis. 


— 
— 


o =e, The sampling apparatus designed and built for this study 


appears to operate satisfactory. The data obtained correlates 
with work done by other researchers and with the computer model. 


However, additional work is needed to correlate sample size and 


particle size distribution with freeboard height. 


Recommendations 


1. Further experimental data is needed at higher velocities 
and higher heights above the bed using the sampling apparatus 
designed in this work. This will provide additional data to 
evaluate the operation of the sampling apparatus and the computer 
model. aAiso, a scale capable of measuring quantities of samples 
less than 9.01 grams and a particle removal system which removes 
the particles uith less error from the sample trap is necessary 


for accurate uork. 





- 


— amme a 


جه 


2. A sampling device which has a shorter trap height should 
be designed. The present sample trap is very direction oriented 
and samples particles traveling only in a narrow range from the 


vertical. This would allow more accurate work to be done in the 


splash. zone where particles are more likely to be traveling in 


directions other than vertical. 


3. The computer model needs to have encorporated in it, a 
statistical distribution model for jet velocity and jet duration. 


This will allow the determination of the effects of varying jet 


-=- 
— a 


— velocity and duration on particle distributions in the freeboard. 


4. A correlation for jet velocity and duration as functions 
of bubble diameter, Uo, and Umf should be determined and included 


in the model. 
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APPENDIX A 


Moment of Inertia Caiculations for Paddles 


The calculation of the total Moment of inertia (1) of the 


"==> --=> ” 88ሮ61 5 -ር838 be broken up into three seperate calculations. First, 


o -ቲክፀ moment of inertia_{Il} of-the-two-aluminum cylinders used to 


mount the paddle arms to the solenoid shaft is calculated. The 


second calculation (I2), accounts for the moment of inertia of the 


paddles themselves, uhich are constructed 
-0 iL. epoxy -laminate-(Fig.) Ar-i).<Finally; the 


the hardwood mounting ends cn the paddles 


The-following equations were used to‏ - ۱ سے 


moments of inertia: 


of a foam, basswood and 
moment of inertia (I3) of 


is calculated. 


calculate the mass 


11 - ع‎ ዝር 
2 (Ai) 
4 2E 
IZ E de له‎ ) + U) (PI ti Z t7) 
12 (AZ) 
e 
12 (A3) 


Table 4.1 is a listing of parameters required for the moment 


calculations. 


po 
(በ 
نی‎ 
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Component 


 Bluminum Cvlinder 


ج A‏ سم 
ھب > ې ۰ 


Density A1 {PI= 
Radius [r] 
Height [H] 


Paddles 
Density Foam [P1] 


Density Bassuood [PZ] 
Thickness Foam [¢1]=.4 


 Thickness-Bassuood [t2] - 


Length Foam (L1] 
Length Basswood [LZ] 
Width Foam [W1] 
Width Basswaod (WZ] 


Hardwood 


Density Maple [P] 
Thickness Maple [t] 
Length Maple [L] 
Width Maple [WU] 


Listing of paddle components and parameters. 


ISI ET 


15: 


MKS units 
2780/— kg/m" 
0.0171 m 
0.0381 m 
56. kg/m 
| 566986 -— kg/m 
0:0102 m. 
8.00153 m 
0.5842 m 
0.5842 m 
0.0305 m 
0.0305 m 
556 kg/m 
0.0102 m 
0.0432 m 
0.0348 m 


English units 


0.1004 
8.675 
IS 


0.090181 


0.0144 
0.4 
0.0625 
23:0 
23-0 
122 
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3 
lbf/in 
in 
in 


E 
ibf/in 


ings 


Hn 


in 
in 
in 
in 


1bf/in 
in 
in 
in 





equations (4.1),-.(A.2), and (A.3), the value for each inertia 


Equation (A.1) is used to calculate the mass moment of 
inertia for a right circular cylinder rotating about its Z axis 
(Fig. A-2). The height H, accounts for both the upper and the 
lower solenoid shaft cylinders. Equations (8.2) and (8.3) 
determine the mass moment of inertia for a rectangular prism 


rotating about its X-Y centrcidal axis: (Fig. A-3). The length L 


in equation: AZ-accounts-for both the upper and lower paddle 


lengths. The factor of Z in equation (A.3) is because both the 
upper and lower hardwood sections are identical and can be 


combined. By inserting the ressective values from Table A.1 into 


component can be calculated. Adding these components, the total 
mass moment of inertia applied to the solenoid by the paddles is 


determined. These results are listed in Table A.Z. 


3 
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۱ 
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Figo? Diagram for Moment of Inertia Calculation 
Used for Cylinder About Z Axis. 





Fre. A-3 Diagram for Moment of Inertia Calculation 
Used for Rectangular Prisim About X Axis. 
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PABLE A.Z 


Moment of inertia components: 


MKS units 


احم 
محم 
| 


= 1.422 (-05) kg m 


i 
N 
ti 


۳ ت2 04 - ۳د 52:92 


Se Ru WS 


English units 
4.911 )-02( in 
1.982 (00) 15 in’ 


6-918 ۶03 1۵ 111 


Total moment: ፲ = ፲1 ተ 12 + 13 


I1 = 5.991 (-04) kg سم‎ 


Calculated values for-mass moment of 


2.037 (88) lbf in 


inertia of 


=-paddles and.paddle.components. Il is inertia 


value for the aluminum cylinders- — 12 is the آ‎ 
inertia value for the paddies themselves and I3 


is for the harduood end pieces. 


Poa 





APPENDIX ኩን 


Error Determination of Vacuum Collection System 


— 


Ul 


acc =Tests-of the vacuum sample removal system indicated that some 
i- i= particles remained in-the sample? trap. These particles were 
located in the corners of the trap where the equalizing air stream 


could not agitate them enough to move them into the vacuum stream. 


As a result, it became necessary to determine to what extent these 


ኤ > 


2: =. su. remaining particles. affected: the: accuracy of the sample attained. - 


The procedure used involved placing samples of known weight 
>> ፦- and particle size distribution (equivalent to the bed material) 
inside the sample trap. Sample sizes of 10, 15, 20, and 25 grams 
were used. The sample particles were then vacuumed out and their 
weight determined. The difference in weights of the samples were 
then calculated along with the percentage differences. These 


results are listed in Table 8.1. 


An analysis of the results listed in Table B.1 indicates that 
the average difference betueen the tuo sample ueights is only 0.52 
x wnile the maximum difference observed was 0.33 X. It can 
sherator2= 5e concluded that the particles remaining in tne sample 
trap do not significantly affect the accuracy of the sample 


attained. 





Initial Sample Weight Weight 
SAMPLE Weight Weight Difference Difference 
No. grams grams grams Percent 
1 19 Sms 0.95 9.50 
2 19 gag 8.98 0.80 
3 45 10 9 95:7 2.65 2.59 
4 19 SIN 0.92 0.70 
5 15 PS. US 20.03 -0.20 
E 15 - 15.02 2 -0.02 20. L3 
7 1: 14. 85 9.14 0.93 
8 5 14.87 0.13 0.87 
9 Dt 14.97 8.03 9.29 
13 2 بعرو‎ 20 =--- 13.9 3. 8.07 0.35 
ہے رنہ‎ 20 s 1.6 | 9.14 9.70 
په بس‎ E ረ ጂጂ چک ور 9077 0:109 19.908 ےا‎ 
١ رت 20 > د‎ 197233 0.97 5 
¿q 29 19.94 0.06 0.30 
15 25 24.83 0_7 6 8 
Iw —— 23 24.80 9.29 0.80 
ld 25 24.78 2522 0.88 
18 5 24.83 8.17 0.68 
18 Z5 24.77 0.73 0.33 
AVG= 0.19 0.5٤ 
STNUV- 9.08 0.34 


Results of vacumn sample removal test. 





APPENDIX C 


Sample Trap Closure Time Test 


Due to the high velocity.of.the particles uith respect to the 


— — — 


size of the sample trap, it was important to ensure that the 
sample trap was closed in a very short period of time. To measure 


the closure time, the set-up shown in Fig. C-i was used. 


A, ናር] 5 ተመቹ ን ርስቲ iti ee ae > 
photo transistor; was used. The diode was powered by a6 volt 
battery and the photo transistor was wired directly to an 
oscilloscope. The electric eye circuit was then placed at a 
distance of 0.267 m (10.5 in) from the pivot point of the paddles. 
At this position, the closing time calculation involves only a 
simple proportion relationship between the paddle and the sample 
trap uidths. The eye uas also positioned as close to the trailing 


edge of the fully closed paddle as possible to minimize error. 


The oscilloscope sueep uas set to trigger off of the 
initial change of state from the photo transistor when the paddle 


nt traces were then 


(v 


first eclipsed the light beam. The resul: 
analysed to determine closure time. Time ቲ=0 was set equal to the 
initiation of the trace on the oscilloscope. Time t-t! uas 


defined as the point at which the trace begins its excursion back 
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ہہ 
Sie‏ 


(a) Too View 


6 0 ۲ Diagram of Closure Time Determination Set Up. 
Ali Dimensions in cm. 


دا 








to the base voltage as shown in Fig. C-2. By using the ratios: 


Wi = WZ 
1:7! t2 CCL) 
Where: 
Wi = width or the paddle 
SSI SANE - uidth ot sample trap_—— —— 


= +232 seis omc ti- time delay measured on oscilloscope 
- ጨም © = tZ = closure time of sample trap 


the closure time of the sample trap can be determined. Several 
trials were run with the resulting data listed in Table C-1. All 


of .the trials were-very:consistantox-The resultant determination - 


for the average closure time was 1.44 ms. 


۱ 1 
Wi = 8.0305 m CIZ TOD 
WZ = 0.0127 m (0-5 in) 
Trial +1 t 
Na. ms ms 
1 A 1+137 
2 524 1.42 
3 3.55 1.48 
4 3.5 1.46 
5 3-5 1.46 
Avg 3207 1.44 


Sample treo closume. data: Wi is the wrath of 
the paddle, WZ is the width of the sample trap 
and time tl is eclipse time of paddle through 
light beam. Time tZ is closure time of Sample 
trag as calculated using equation (C.1). 
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Time (8.2 ms/Div) 


Oscilloscope Trace of Paddle 


Eclipsing Electric Eye. 


Gaz 
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APPENDIX D 


Solenoid Torsue and Dynamic Analvsis 


= „Ihe pond of this design uas. ‚sag shut the sample trap in about 


د un‏ ہت ےس m‏ رت موند 


l millisecond (ms). At this speed, a particle traveling at 10 
meters per second (m/s) would travel 1 cm or approximately 11 % of 


the sample trap length. A velocity of 10 m/s is at the upper 


= limit of the particle تتت س‎ and would account for 


— — ست‎ = መ፦-መመ-ሙ፦-- - ጨ።ሙውው-- -----መ 


B چیا‎ 


a small percentage of the particles ee by the JR * The‏ حا 


- majority of the particles, for the fluidization conditions used in 


the bed, have an average ejection velocity of 1-2 m/s, based on 


data from George and Grace [81]. 


TORQUE ANALYSIS 


To determine the minimum required torque output of the 
solenoid, two analysis were performed. The first analysis 
calculates an average torque required to produce the desired 
velocity. The second analysis uses the torque data from the 
selected solenoid and calculates the expected closure time and 


swing time for the paddle arms. 
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First Analysis 


To close the sample trap in 1 ms, the anqular velocity of the 


paddles at closure can be calculated as: 


$ وا‎ 2 ህረ 
سی‎ ٢ 2 کات و سس‎ 
where: w = Angular velocity of paddles 
جج :ددد کے‎ sais ها ند‎ = Width of sample trap. 
E = Ss_sarRige Paddle:pivot: to sample trap distance. 
m Ana کش‎ --t2-2-Closure time of sample trag - E 


(D.1) 


Using the values of W2= 0.00187 m, R- 0.0413 m, and t2= 1 ms, 


equation (D.1) results in an angular velocity of 47.70 rad/sec. 


The radial acceleration required to attain this velocity through a 


deflection of 4 radians (45 Deg), is calculated by: 


= 2 
a = وا‎ 
2 0 
wnere: a = Angular acceleration of paddles 
w = Angular velocity of paddles 
8 = Angular deflection of paddles 


GO) 


Using the result of equation Dl and 823574 radians, equation DZ 


gives the required angular acceleration of the paddles as 1448.5 
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rad/s/s. Combining this result with the result for the total mass 
moment of inertia from Appendix A, the required average torque 


output of the solenoid can be calculated as: 


ہے سے 3 >= ው.‏ 


Met سو "ې سو‎ = Angular acceleration of paddles | 
Mass moment of inertia -يَ1_ٹ-ٌ 0 یب‎ 
= Required solenoid torque 


= 
I‏ = = یس سس کے ee‏ ل کار A‏ 
T‏ 


— fos t - በጨ. 


The resultant value for the torgue (T) is 0.87 N m (7.7 Lbf in). 

î A geometric „average of the. rat soienoid' s torque output, as 
-+= -—-shoun in Fig. D-1a; indicatesTan average of 1.8 Nm (8.8 Lbf in). 
Therefore, the selected solenoid has sufficient torque output to 


achieve the desired paddle velocity. 


SECOND ANALYSIS 


For the second analysis, it will be assumed that the torque 
output of the selected solenoid can be modeled as a linear spring. 
Fig. O-lb shows the torque output of the solenoid as a function 
of angular displacement. From Fig. 0-1۳, a spring constant of k= 
B 0.72 N w/rad (0.11 Lb5£. in/deg).can be used to approximate the 
torque curve. The spring constant is derived from the slope of 


the iorgGue-angular displacement curve. 
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The system can then be modeled as a simple rotational mass - 
spring system. The differential equation of motion, general 


solution, and boundary conditions are: 


qi 
+ 
- 
ጩ 
ዘ 
ሠ 


Df Egn: I 


(0 da) 
E= "Gen Soln= 0 -8 2-C1 sin(m t)' * C2 cos(m t) "CU. 4D) 
-— +-ے۔ ہے‎ — where: m = k u NR بجع‎ ER E NI سح سے ا ل‎ m 
1 
ےن‎ 1) + =: 8-0 
2) t = 9- 8 = 1.65 rad 


۔ جم 


z- Solving equation (0.4b) using the given boundary equations, gives - 
-the follouing expressions for any angular displacement (80) and 


angular velocity (8): 


8 = 1.65 cos(m t) (D.5a) 


-1.B5 m sin(m t) (D.€b) 


me 
il 


Knowing that the solenoid stroke is 7/4 radians, the total solenoid 
actuation time cen be obtained from equation (D.Sa) by setting Û = 
(1.65 - 7/4). The total actuation time was determined to be 29.2 
ms. Inserting this value for time into equation (0.Sb), a 
rotational velocity or 48.5 rad/s was found. This velocity is 


very close to the desired velocity obtained from equation (0.1). 
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DYNAMIC ANALYSIS 


The desire to determine whether the debris of a specific 
bubble has been trapped by the sampling device requires the 
determination of the total actuation time for the apparatus. 


- 352: nouing. this, value, . the delay betueen- the detection of a bubble 


~ 


E and the closure of the sample trap can be determined, and hence, 


_whether_the_particles from_the_detected bubble were within the 


vicinity of the trap. 


$a === Appendix. determined that the actual closure time-for the 


~ 7 
— کب ما‎ — we = 
- 3 -r-e መሙ 


hic al En equation ን ዜር‏ حم حت 1.44 sample trap was:‏ چم 
=z +Erotational velocity at the end of -the ጆሮ of only መ te ace‏ 
-is calculated. Because the mass moment of inertia for the‏ 
solenoid was unknown in the beginning, the ቸው ignored it.‏ 
With the determintion of the actual closure time in Appendix C,‏ 
the total moment of inertia, including the solenoid can be‏ 
approximated. By iterating equations (0.Sa) and (0.56) with 8 =‏ 
radians/sec, a neu inertia value of 1.27 (-@3) kg m is‏ 33.12 
determined. Using this inertia value in equation (D.Sb), the‏ 
total actuation time is 42.6 ms. This increase in actuation time‏ 
of 13.5 ms is relativily close to the manufactures quoted‏ 


actuation time for the solenoid of 12 ms and suggests that the 


values are within reason. 


Table 0.1 provides a summary of the results obtained from the 


dynamic analysis. 
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op 1 


Mass moment of inertia: 


Dee 5.991 (-04) kg n- 
2.0537 (00) 1b in“ 


= o Paddles + Solenoid. (It) 1.27°(-03) ka m- 
۰۳-۳ ااا جن‎ 4.32 (00) lb in” 


— 


7 » 
m =‏ © = = > ہے ہے ې 


Closure time: 


Sample trap (t2) 1.44 ms 
Total suing (t3) 42.5 ms 
جرته خم‎ To Summary of dynamic analysis results. 
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ATTEMPT E 


Bed Particle Size Distribution Analvsis 





A sample of particles was obtained from the fluidized bed 
after the bed had been operating for several hours. Sample sizes, 
ranging from 1.2 to 2.8 kg, were then taken from the central area 
of the bed. Each sample was then sieved through a series of 14 US 
standard wire mesh sieves using a Tyler Industrial Products Model 
RX-24 portable sieve shaker for Z0 minuites. The contents of each 
اب‎ was then weighed, using a Torsion Balance Co TORBAL scale, 
to an accuracy of 0.01 gram. The resulting average particle size 


distribution is listed in Table E.1 and shown in fig. E-l. 
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TABL 


ha 


m 
m 


US Seive Ueignt Sercent Cumulative 
Standerd Size In Seive ot Total Percent 
Seive No. um grams Sample 

18 1898 8.49 071 01 

290 aca 24.59 3.59 121 

26 630 181.26 sld 4.35 

3 500 119.85 4.06 3.01 

45 3S3 176.15 6.46 1S.46 

ca 257 430.77 15.78 31.24 

23 250 565.54 14.12 45.5 

.7178 Zig 6 2 12.48 == . 84 
809 139 453.29 17.18 77.03 
190 149 ووي‎ 83 12.45 89.48 
120 ES 125223] 4.58 ° 24.06 
149 196 198.45 4.91 38.25 
179 3g 46.57 ٢۸ gg 7 
200 75 5.91 Ooze 997958 

(7S 9.38 9.91 108.98 


Particle Specific Gravity = 8.1 


Average particle size distribution of bed material in grams 
and percentage of total weight using sieve analysis. 
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F‏ گر و یہ جح بب تہج 


Mean Bed Flow Velocity Determination 





To determine the mean velocity of air flouing through the 
bed, the ASME report on fluid meters [17] was used. Based on this 
paper, equation (F.1) was used for determining the mass flou rate 


of air through the bed. 


= E 
لا‎ 2 0033702 ) እ ۳ 6 "(e dP ) MEE 


where: 


= Mass flow rate (lbm/sec) 

= Orfice diameter (inches) 

= Flow coefficient 

Thermal expansion factor (1 for air) 
= Expansion factor 

= Differential pressure (inches water) 
= Density of air (ahead of orfice) 


۶ 71 -> 7۲۱ 2 € 
1 


The flou coefficient (K) is a function of Reynold's number, the 
orfice diameter (d), and the pipe diameter (D). For the pressure 
tap configuration used in the MIT atmospheric fluidized bed (1-0 / 


1/2-D) the flow coefficient is calculated using equation (F.2). 


Hæ 





where: 


Re 


Ko 


0.000 5 FB 


Ko + 1000 b 


VE Re! (ዮ.ረ) 


Reynolds No. 
d/D = (orfice dia / pipe dia) 


~- 
(9.6014 - 0.91352 D /4 ر‎ 


a, 
(8.3768 + 0.07257 0 © 


e 
A 


o? 8^«0.00025 D 


(0.0007 + 9.9011 ) 
D 


(0.0038 + 9.9994) 
D 


^ 


z Je 
A 68 ) 


The expansion factor (Y) is determined from equation (F.3) and is 


a function of the diameter ratio (B), the ratio of specific heats 


(s), and the ratio of the differential pressure to inlet pressure. 


H 


Y= 1 - (0.419 + 0.350 B ( 8ل‎ 


where: 


Pls ا‎ 


Pi Inlet pressure 


Ratio of specific heats (1.4 for air) 
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Mass flow rate is a function of velocity, and therefore, by 
definition a function of Reynold’s number. Equation (F.4) is used 
to determine mass flow rate as a function of viscosity and 


Reynold's number. 





C Re (F.4) 


where : 


u Z Absolute viscosity (lbm/ft sec) 
Constant (15.28) 


e 
u 


This requires that the solutions of equation (F.1) and equation 
(F.4) be iterated until the Reynold'’s numbers converge. The 
determination of mean air velocity is obtained from the mass flow 


rate (W) using equation (F.5). 





cun CES) 


where: 


3 
= Density of air در‎ ) 
Area of bed )++ ( 


=> ۵ 
1 | 


1:53 





To perform these calculations, a computer program was used. The 
program (APPENDIX G) is written in HP BASIC 2.0 foc running on an 
HP 39816 series 200 micro-computer. Convergence usually required 


four to seven iterations and should be correct to within 0.5X. 
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APPENDIX G 


This Appendix iists the BASIC computer program used to 
calculate mean bed velocity as discussed in Appendix F. It is 


written in HP BASIC 2.0. 
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10 
28 
20 
20 
50 
50 
70 
82 
30 
106 
110 
sr atm, 


^n 
4 


138 
123 
189 
159 
179 
189 
138 


pp) FP) ss 
کی مس ل‎ 
O © © 


አህ ዮ)ኔ f‏ 63 دا دا 
NM YY bi‏ 
نم C ($9 & C‏ 


را 
c»‏ 


is 
10 
© 


C €») (Y نی دج‎ €» 


man 4۶۱ ۱۱ ۷۱۱ الا‎ 
wor ات‎ )(( Ul ۰ 


HA ex MUNDOS TOCA Mc 
Mezerıng of gases Dy means of the ASME souare-ecees orfice with 
{0 1/2-0 taps. Reference Fluis Meters Their Theory and 
25581۱:۳۵ 160۳۰ ASME Repor: bin ecition, ASME New York, NY, 1371 


١ 

! 

۱ 

! Program mus: be altered for Orfice diameter other than 7.071 in 
! and pipe Giametec otner than 10.02 in. (see line !30) 

۱ 5 

| 

INTEZI2 Ansue” 

ER ዴ: ን. ی‎ 0021. 29652 መመመ 5ed.Umut,Reynolcs no,Rhc,Rhol 
) ሮመ። ላ 

DIM Homesil1,CcCl1ears(2] 

DATA 7.371,18.07,1.4,11.211,50000.,76. 

READ D2,01,S,Aree_2eg,Reynmolos_no,Patim 

Ciears=CHRS( 225 )&8CHRS( 7S ) | CLESR screen 


^cz 


ہے هد ጠ-‏ 


HomessCHRS( 
| 
I Input 
| 


۰ 


INPUT 
INPUT 


INPUT 
! 


JECHRS(BL ) 
vantabl=s 
و‎ E pressure PI (capo): c Pt 


“Enter pressure oron GP? (incnes uater): 
“Emier air temperature (Degrees ዮ)፥ ",Tl 


gata 


I Comaute 
| 

QUTPUT 2;Homes; 

OUTPUT 2; Cheacs: 

Beta=02/01 

Pla=?1+?atm 

Ple(Pla)/Z.34 

7۔ 2ت 2200186 

Rho=rNDensity(Patm,T2) 
Rhol="NDensi:y(Pla,T2) 
Umus-NUiscosity( 12) 
Yl*FNZzxpan(S3et:a,?1,5$,9 del) 
K*FN-lou coeff(Reynolds no.3sta,Dl1) 
ع لیا‎ 2۱۱۳۵ 55 flou(D2,X,Y1,Rhol,? del) 
8۵۶۶۱2۵2 لاب‎ Umu,D2) 

i= 9251 13518۳۳۰155 0۰ 


PZN 


I HOME screen 


Home display 
Clear display 
Celculate Beta 
Absolute pressure 
Change to inches 
Convert to R 

Air density bed 

Air density upstream 
Bic viscosity 
Expansion factor Y 
Flow coefficient K 
Mass flow rate 
Reynolds 2 calculated 
Check accurac 


ነጎ “”መ در سار‎ == ጎ፥+” oan r ለ EEE gem መሙ ፥መሙ፣  ፥መው یر‎ 


Revnolds_no=Rd 
5070 393 
END IF 
Veloc::ty=W/(Rho+Area_bed) ! Yelozity ceiculation 
PRINT "Veiocwiy= ";Velocıty:" ft/s° I Qutaut results 
INPUT “Enter (1) to continue, (@) to stoo” ‚Answer 
TF Answer=@ THEN ۳ 
Ir Ansuer«»1 THEN 
SINT "Value must 3e a:ither 1 or 985 try ۵ 
S010 0 
END 2 
GOTO 6 
ENO 
۱ 
| 2121217 کم همهم‎ S2gin FUNCTIONS 
l 
m No camt( seta. > t.2.7_ae!l) 5 م کے‎ TY 
ንኛ ር رط 2 2 )/ صا‎ 
۷۱-2 35866 (0 9 
RETURN Y 5 





س نی 


630 FNENO 

siQ [ 

Se 

226 DEF FNFiow_coe?f(Reynmolds_nmo,Beta,Dl) [Flow coarficient K 
549 Koz.8014-.201222*01^(-.25)*(.5378*.272587*01^(-.2982)0«(.290025/(01*01*8&e:a*2e 
+ ۰و ۶ 5:۰2 ۳41 32+( 25۰101 ۵00 جج‎ 15 ( 

520 2=.2002+.30!::/D1-( .3033+.3004/091 )+(Beta-Betar(16.5-S-D01 )#Betia if) 
2529 ۶۱0+ ۱000۰53/52 ) 565255۳0 1 5 0 ( 

270 SETURN K 

2239 =NEND 

522 [ 

753 [ 

713 BES ساود بيرع‎ የ!16ህ(:2,ሮ.ሃ1, 51601,ዎ 251) ! Mass flow rate Y 

(20 Ux.292;702*02-*02*K«Y1*SQR(ORRnoOl1*P del) 

730 RETURN Y 

720 FNEND 

7-0 [ : 

760 [ 

772 DE? FNRayn(U,Umu,D2) ! Reynolds + calculatad 
182 32 0۶۱ 259ሀህ/ ('ህጠህ*ቀ02) 

759 RETURN 1d 

20 FNEND 

319 l 

820 [ 

8:39 DE- FNOensity(Pstm,i2) I Calculate air density 
840 Rhoz.z222405*Paim/TZ 

253 RETURN Rho 

9859 FNEND + 

378 [ 

380 [ 

226 DES FNViscositv(TZ) l Calcaglata àir viscosity 
220 Umus7.52058E-7*12^1.2/(T1T22198.5) 

219 RETURN Umu 

329 FNENO 





۴5۱۳۲ ۳ ۱۵ 1 ንጋሩፈ H 


This Appendix contains the parts list for the sampling system 
in three tables. Table H.1 is a list of components for the sampling 
apparatus. Table H.2 is a list of components for the electrical 


Circuit while Table H.3 is for the vacuum system. 


178 





item No. 


Component 


Solenoid 


Paddle 


Interfaces 


Solenoid 
Mounting 


Extension 
Bar 


Trap 
Mount 


Sample 
Trap 


Vertical 
Mounting 
Slide 


Base 
Structure 


Table H.1 


Description 


LEDEX Size 5S 45 Degree Right 
Hand Stroke Rotary Solenoid 
Part No. 5-8204 -9 

LEBEX Inc. 

801 Scholz Or, P.O. Box 7 
Vandalia, Oh 45377 

(513) 898-3621 


Foam, Basswood epoxy laminate 
with Maple mounting blocks 
0.59 X 1.2 X 11.5 inches 


Aluminum Cylinders uith set screu 
fastener. R= 9.675 in H= 1.5 in 


Aluminum mounting plate 
1/4 X 2.8 X 2.8 inches 


Aluminum bar 
1/4 X 3/4 X 12.6 inches 


Aluminum Plate 
1/4 X 1.1 X 1.5 inches 


1/16 inch Aluminum plate 
Inside Dimensions: 
3.5 ደ 1.5 X 3.5 inches 


Aluminum bar 

1/4 X 3/4 X 24 inches 
Adjustment holes drilled every 
0.5 ٨ 


Tripod Aluminum Structure 


List of components for sampling accereiíus. 





Item No. Component 
1 91 
2 Di 
3 M1 
4 R1 
S R2 
5 R3 
1 51 
8 52 


TABLE H.2 


Description 
Bridge rectifier assembly, Silicon, 
LEDEX Part No. 121011-001, includes 
arc suppressor in unit. 
Arc Suppressor, not needed if above 
rectifier assembly used. 
LEDEX Part Ne. 122555-01 
Solenoid (see Table H.1) 
Resistor, Z Mohm, 1%, 2 Watt 
Resistor, 100 kohm, 1%, 1 Watt 
Resistor, 250 ohm, 10%, $8 Watt 


Switch, SPST, Push Button, 
10 8, 250 V 


Switch, SPST, Toggle, 10 8, 250 V 


List of components for solenoid power supply. 
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as MEN ٣٢ 


00, 0۲ 


5 


E. 


5 





Item No. 


TABLE ۵ 


Component 
vi 
P1 
F1 


ET 


€Z 


ET 


Description 
3/8 inch Ball Valve 
Eductor 
Screen Filter, 320 um mesh 


Sample Container, Small Plastic 
86 ماع‎ D 


Sample Trap (see Table H.1) 


1/4 inch Polyflow Tubing 


List of components for vacuum system. 
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APPENDIX I 


This Appendix contains a complete listing of all data 
obtained during this study. The data is arranged 
acording to the distance of the trap heignt above the 
distributor. 


q 
ቕን 
rJ 





DATA for Umr DETERMINATION 


Hb P1 Pb ር ፐ 
(inches) (cm Hg) (cm H20) (n H20? (F) 
18.75 T 51.8 1.805 62.0 
19.75 5.8 50.9 1 7 62.0 
19.28 6.2 49.1 0:937 62.0 
3.75 5.8 43.8 8.772 52.86 
3.59 5.8- 47.5 0.515 62.0 
9525 5. 1 44.8 0.335 62.5 
8. 5 b 4.8 42.83 8.265 63.0 
8.59 4.5 |. 89.8 08.1 77 63.9 
8.59 4.0 29.5 8.125 — 63.5 
8.50 3.7 33.9 9.197 . 54 . 0 
8.50 -. 5. 95 خم‎ 32.0- 0 . 4 54.0 
8.58 3.8 30.1 90.084 64.0 
8.50 fo m .4 9.943 64.5 
8.52 2.9 18.5 0.027 65.9 
8.59 1.5 


14.6 0.317 609 
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Trap Height above Distributor: cm (in) 50 36 (12.13) 


Bed Height above Distributor: cm (in) = 25.4 (19) 
Trap Height above bed surface: cm (in) = 5,56 (2209) 
Pl: cm Hg (in Hg) = 5.8 (2:28) 
dP: cm water (in water) = 1.234 (0.486) 
Temp: C (F) = 19.5 (67.9) 


SAMP 5 Weight of Sample 
(grams) 


20065 
.-- 2.06 
1288 
2.43 z 
2.48 
2.10 
1.84 
2 0% 
2-7 
16 2.15 
11 ء2‎ 3 
2 1.83 
15 2.46 
14 1.13 AVG 
15 1. 1 STN DV 


مرا رم ين ۵ 0 ن0 [- ين ف 


‚2.893 grams 
9.50 grams 


38-98 112213) 
25.67 (1025) 
4.29 (1169) 
6.8 (2.868) 
1.887 (9.743) 


Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
Pl: cm Hg (in Hg) 

dP: cm water (in water) 


Temp: C (F) 18 (64.5) 
SAMP i Weight of Sample 
(grams?) 
16 2 
1717 93992 
18 .د‎ 3 
9 4.94 
20 5.5 
21 3.03 
22 8.86 
23 4.50 
24 5 5 AVG = 5.96 grams 
25 5.84 STN OV = 2.12 grams 
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Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 


Pl: cm Hg (in Ha) 
GP: cm water (in water) 
Temp: C (F) 


- SAMP- E Weight of Sample 
(crams) 
25 0.82 
27 2035 
2 28 9.54 
IA 29 9.54 
26 9.76 
31 9.58 
32 9.51 
33 90.47 
24 0.85 AVG 
35 9.51 5 [ከ 
Trap Height above Distributor: cm (in) 


Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 


PY: cm Hg (in Ho? 
GP: cm water (in water) 
Temp: C (F) 
SAMP # Weight of Sample 
(grams) 
25 2-87 
37 1.42 
38 1.89 
55 1.34 
46 ۱ ول‎ 
41 1 2 
47 ۱۱ 
43 1.53 
A4 1.56 AVG 
45 123% STN DV 


13S 


DV 


36.31 -)14.( 
25.4 (19) 
1١ 91 (27853) 
5.9 (2228) 
1.291 (9.473) 
20.5 (83.9) 
= 0.57 grams 
= 0.11 5 
37.31 (14.59) 
26.67 (10.5) 
10.64 (4.19) 
5 ژ‎ (2.28) 
1.897 (0.747) 
183 (65.9) 
1.66 grams 


0.46 grams 





Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
۲۶۱ sor Ho (in Ho) 

dP: cm water (in water) 

Temp: C (F) 


SAMP zx Weight of Sample 
(grams) 
46 4.32 
47 E 2.48 
48 3.12 
£o 49 283 
EMI 28 4.46 
51 237 
57 4.46 
= 53 4 . 5 

54 4.15 AVG 

55 DOSE STN DV 


186 


30. «14:537 
27.94 (11.0) 
gest (3:997 
5.8 (2.58) 
1.887 (0.743) 
i (55.9) 


3.61 grams 
0.95 grams 





TRAP HEIGHT = 29.25 in 


Heignt of Pressure Pressure Pressure Pressure 
Pressure for for for for 
Tap above Samples Samples Samples Samples 
Distributor 66-75 86-95 56-85 76-85 
(inches) (cm HZ0) (cm H20) (cm H20) (cm H20) 
1.8 259 54.9 56.1 59 
2.8 48.8 48.8 SIT 52.8 
3.5 40.5 —— 41.3 43.8 44.7 
4.5 34.93 55 8 SES 40.0 
5.6 27.8 23 7 372-5 .د233‎ 5 
6.8 ረ ጋ 22.53 Ze 5 27.6 
7.6 14.6 15.1 1 20.7 
8.6 f 8.3 Sb d عا‎ 5 14.8 
3.6 ፦ ር ረ. 3 5.5 e 5.4 8.5 
10.5 de? 93.6 3 9.1 = e 2.4 
د رر‎ 9.9 0.0 9.9 0.3 
25 0 . 0 0 . 0 00 3.8 
Les 8.9 9.0 . 00 00 
i4.6 9.0 9.9 9.6 9.0 
TRAP HEIGHT = 23.25 in 
Height of Pressure Pressure Pressure Pressure 
Pressure for for for for 
Tap above Samples Samples Samples Samples 
Distributor 286-215 1236-205 186-135 175-115 
(inches) (cm HZO) (cm HZ0) (cm H20) (cm HZO) 
1.5 52.0 53.4 56.3 7۔57‎ 
2.6 48.9 49.6 51 8 5 
3.6 49.7 4215 43.6 44.7 
4.5 34.7 35.8 38.8 40.9 
5.5 25.3 29.3 32.0 33353 
5.6 21.5 22.9 25.6 25.5 
7.6 14.8 16.5 928 29.3 
8.6 8.4 gB 82 138 
9.5 2-2 2 3 5.2 25 
19.8 2.9 9.1 | 7 2 
| 3 00 3.3 0.92 9.3 
12 2.9 8.2 0 . 0 0 . 0 
15 0 . 0 29.0 2.6 0 0 
14.5 9.9 9.2 3.8 9.0 
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| 4. 


I 


Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
Static Bed Height: cm (in) 

Pi: cm Hg (in Hg) 

dP: cm water (in water) 

Temp: C (F) 


SAMP E Weignt of Sample 
, - (grams) 
76 1.47 
نے‎ 0 ۰٠ 1-35 
78 1.24 
ت یایۓ‎ 09 15 
=> 80 8.74 
۱ 81 1.38 
82 1:38 
83 hn m7 

84 1.58 AVG 

8S 10033 STN DV 


Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
Static Bed Height: cm (in) 

PT: ¿nm Hg (in Hg) 

GP: cm water (in water) 

Temp: C (F) 


SAMP # Weight of Sample 
( grams) 
S6 0.95 
57 0.61 
Ed 15 
59 0 . 4 
59 1.02 
Bi 0.76 
52 0.57 
53 1.95 
64 2.60 AVG 
65 0.88 STN DV 
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51.44 (20.25) 
723 CLER) 
22.23 (8270) 
22.56 (8788) 
323 CID 
56.385 (7.75) 
155 (5220) 


1.32 grams 
0.23 grams 


51.44 (20.25) 
27,31. 19.79) 
22.13, 03.5) 
27.56 (8188) 
ሆን (229m) 
453707 11.333) 
15.5 (62.0) 


0.92 grams 
0.36 grams 





Trap Height above Distributor: cm (in) 


Bed Height above Distributor: 
Trap Height above bed surface: 
Static Bed Height: cm (in) 


P1: cm Hg 
ርም: 


Temp: 


cm water 
ር (ዮ) 


(in Hg) 


SAMP 5 


86 
87 
38 
8g 
50 
0 
32 
33 
34 
35 


(in water) 


cm (in) 
cm (in) 


Weignt of Sample 


Trap Height above Distributor: 
Bed Height above Distributor: 
Trap Height above bed surface: 
Static Bed Heignt: 


P1: 
dP: 
Temp: 


cm Ho 


cm water 
C (F) 


(in Hg) 


SAMP # 


66 
67 
68 
E3 
73 
21 
72 
3 
74 
75 


(grams) 
6.67 
9.08 
8.17 
9. 4 
0 . 8 
۵ . 4 
0.97 
013 
7۔0‎ AVG 
om STN DV 
cm (in) 
come Cin) 
cm (in) 


cm (in) 


(in water) 


Weight of Sample 
(grams?) 


0 . 5 
0.5 
0.14 
2.07 
GT 
0.10 
9.18 
0 . 8 
0-11 
0.16 


18 


AVG 


STN DV 


3 


51.44 
26.57 
24.76 
272.56 
5-5 

3١. 1 
16 


(20:295) 
(10.5) 
(3275) 
(8.88) 
(2.55? 
(1.284) 
(61.0) 


0.09 grams 
0.05 grams 


51.44 
25.4 
26.03 
22236 
6.2 
2. 1 
17 


0.12 


(20.25) 
(10.0) 
(10.23) 
(8.88) 
(2.44) 
(1.018) 
(63.0) 


grams 


0.04 grams 





Height of 
Pressure 
Tap above 
Distributor 
(inches) 


b org معو هدخ‎ E 


Height of 
Pressure 
Tap above 
Distributor 
(inches) 


مح صو 
ኮኪን re‏ 


موا مو 
Ul‏ £ 


UNE‏ خ O {N‏ ل- 00 X(‏ وخ مع IN‏ )0 ھ 


፦፦ 
GQ OO -J) 0( U1 & مر یم یں‎ 


OOO CU OCG OO OO O CO O) 0) 


OO O CI COO OO OOO O O ጩ 


TRAP HEIGHT = 


Pressure 
for 
Samplas 
186-115 
(cm H20) 


፦ FJ nuu ظ‎ a sy 
SY Y 809 on Ww N 01 لف‎ Y) 


لم 60 ጩ -1ቱ ጩ‏ 1( لد م۳ يہ ہت تہ مت تہ 


Press 

for 
Samol 
125-1 
(cm H 


፦ ዮት) ዮ) ( 5 A 0) 


ظ ص — 1( OA PD (to‏ ما ہیں و و و هی ہت 


TRAP HEIGHT = 


Pressure 
for 
Samples 
146-155 
(cm H20) 


5 ۰ 
49.8 
41.2 
34. 
28. 
Zr. 
14. 

ን 


S SO SO C» C») e Ul to Ul © O) 


Eo oo 6 PJ) 


C C YO بن بت‎ F9) & UN 0 mida 00 (0 fF 


15.12 in 
ure 
as 


55 
207 


18-26 


Pressure 


tor 
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15 
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i‏ 
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1 


- 


Pressure 
tor 
Samples 
155-185 
(cm H20) 


55.3 
1 9 
42.3 
38. 
31 
24. 
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Pressure 
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Samples 
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(cm H20) 
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Pressure 
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Samples 
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(cm H20) 
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44 9 
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33.3 
ZEE 
20. 
14. 
T 


G C ته‎ Ui Ul! 00 ወጠ] ሀ) 
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Trap Heıght above Distributor: 
Bed Height above Distributor: 


cm (in) 
cm (in) 


Trap Height above bed surface: cm (in) 


Static Bed Height: 
P1: cm Hg (in Hg) 
dP: cm water 
Teme: C (F) 


116 

1٢ 

7 118 
5.57 TE 
. ۱ 128 
121 
122 
123 
124 
25 


Trap Height above Distributor: 
Bed Height above Distributor: 
Trap Heignt above bed surface: 


Static Bed Height: 


PP: cm Hg (Cin Hg) 
GP: cm water 
Teme: C (F) 

SAMP + 


36 
57 
58 
58 
106 
11 
192 
103 
184 
15 


cm (in) 


(in water) 


Weight of Sample 

(grams) - 
. 8 
. 3 
2271 
.47 
.34 
. 54 
73 
۔‎ 2 
. 6 
. 38 


AVG 
STN DV 


FJ — rJ r- rJ مار يرم‎ FJ مھ بن‎ 


ca (in) 
cm (in) 
cm (in) 
cm (in) 


(in water) 


Weight of Sample 
(grams) 


15: 39 
9.81 
| 37 
SB 
1.42 
1.48 
1.7 
2.55 
۳ 
0,1 


AVG 
STN DV 


121 


40.35 
28. 5B 
20 37 
L256 
8.6 

5. 858 
15 


2.41 


015.12) 
0112,25) 
(4.88) 
(8.88) 
(Sale) 
(ረ... 
(61.9) 


gr ams 


0.76 grams 


40.35 
227.30 
13265 
22.56 
7.4 
4.808 
16 


15.12) 
(10.75) 
( 5,9989 
(8.88) 
(ረ 1) 
(12833) 
(61.0) 


1.36 grams 


22327 


grams 





cm (in) 
cm (in) 
cm (in) 


Trap Height above Distributor: 
Bed Height above Distributor: 
Trap Height above bed surface: 
Static Bed Height: cm (in) 


PI: cm He (in Hg) 
GP: cm water (in water) 
Temp: C (F) 
. po SAMP. E _ - Weight. of Samole 
۳ ۳ E 0 (grams) 
. بت‎ .- --~126 0 . 5 
7 127 0 . 1 
E21 128 221 
EU |) ን 
126 0.18 
131 e 
132 0.354 
133 0.29 
134 0.34 AVG 
135 9.15 STN DV 


Trap Height above Distributor: 
Bed Height above Distributor: 
Trap Height above bed surface: 


cm (in) 


cm (in) 
cm (in? 


Static Bed Height: cm (in) 


Pl: cm Hg (in Ro) 
dP: cm water (in water) 
Temp: C (F) 
SAMP # Weight of Sample 
(grams?) 
106 3.38 
197 0.12 
138 3.15 
9 0.17 
113 9.13 
lli 4.18 
112 3.58 
113 0.05 
114 0.18 AVG 
13 0.24 STN DV 


2 


(LBZ) 
(1935) 
(5362) 
(8.88) 
(2.50) 
(1.283) 
(51.0) 


0.28 grams 


0.08 grams 

40 95 (16.12) 
26.04 (19.25) 
14.91 (5.88) 
22.56 (8.88) 
6 3 (2.48) 
2.591 (1.020) 
16 (61.9) 
@.22 grams 


0.15 grams 





Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 

Trap Height above bed surface: cm (in) 
tatic Bed Height: cm (in) 

Pi: cm Hg tin Hg) 

CP: cm water (in water) 


46506, (18-25) 
28.58 (11275) 
17.78 (7.00) 
22.56 (8.88) 
8.1 (3.19) 
6.358 ٤ 


Temp: ር (ዮ) 18 (61.9) 
—-——- SAMP. $ Weight of Sample 
De ۳ بصع‎ = (grams) kD 
166 1213 
167 - - 1.18 ቲ > 
158 1. 34 
169 2.38 
179 9.88 
P 171 1.21 

172 1.20 
173 1.28 
174 1 AVG - 1.08 grams 
175 9.85 STN DV = 2.16 grams 


46.36 (18.25) 
27.30 (10.75) 
13.06 (7.5) 
22.56 (8.88) 
43 (2.87) 
4.808 (1.893) 
16 (61.0) 


Trap Height above Distributor: cm (in) 

Bed Height above Distributor: cm (in) 

Trap Height above bed surface: cm (in) 
static Bed Height: cm (in) 

Pl: cm Hg (in Hg) 

dP: cm water (in water) 

Temp: C (F) 


SAMP # Weight of Sample 
(grams) 
156 0.64 
157 0.48 
158 9.78 
158 0.58 
159 9.44 
ioe 0. 9 
162 9.83 
163 0 57 
164 A252 AVG = 9.68 grams 
"BS 0.62 STN DV = 9.12 grams 








Trap Height above Distributor: 
Bed Height above Distributor: 
Trap Height above bed surface: 


cm (in) 
cm (in) 
cm (in) 


Static Bed Height: cm (in) 


P1: cm Hg (in Hg) 
GP: cm water (in water) 
Temp: C (F) 
o SAMP t Weight of Sample 
ን rane) آ‎ 
136 - - 0.03 
157 0.13 
158 0.26 
139 0. 5 
149 9.97 
141 0.06 
142 0.85 
143 2. 65 
144 0.11 AVG 
14S 2.14 STN DV 
Trap Height above Distributor: cm (in) 
Bed Height above Distrıbutor: cm (in) 
Trap Height above bed surface: cm (in) 


Static Bed Height: 


PT: em Hg (in Hg» 
dP: cm water 
Temp: C (F) 

SAMP # 


146 
147 
148 
123 
150 
1S1 
2 
13 
154 
5 


cm (in) 


(in water) 


Weight of Sample 
(grams) 


0 . 1 
0.01 
0 . 2 
2.02 
0.03 
0.22 
0 . 1 
0 . 1 
0 . 2 
0.05 


AVG 
STN DV 


194 


ዚ ዘ ዘ ዘ ዘ | اا‎ 


46.356 (18725) 
25657 (1939 
15 ۹3 ro) 
2256 (8788) 
6.5 (2.60) 
3.244 (1.277) 
1523 (82.8) 
0.08 grams 


0.94 grams 


46.36 
26.04 
20.32 
22.656 
5.3 
2.1 
(65 


18223) 
19225} 
(8.88) 
(8.88) 
(2.48) 
012038) 
(62.8) 


0.02 grams 


9.01 


grams 





ሺጺ..ፕ ፌ-ኤጨሙ ومو‎ 


Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trao Height above bed surfaces: cm (in) 


Static Bed Height: cm (in) 


P1: em Ha (in Hg) 
GP: cm water (in water) 
Temp: C (F) 
٤٥ Same 2 Weight of Sample 
" (grams) — 
176 0-32 
17% ቪር 
178 9.24 p 
178 0.52 
189 0. 22 
181 0.42 
182 9.34 
183 8.31 
184 6). 222 AVG 
185 0.21 STN DV 


Trap Height above Distributor: cm (in) 
Bed Heignt above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
Static Bed Height: cm (in) 


Ph: gm Hg (in Ho) 
GP: cm water (in water) 
Temp: C (F) 
SAMP # Weight of Sample 
(orams) 
. 186 0.25 
187 0.17 
188 0.14 
189 0.17 
129 0.12 
131 0-23 
132 0.15 
133 2725 
194 8.18 AVG 
195 0.17 STN DV 


35 


il O | ا‎ g ክ dH 


50.096 (23.25) 

27.84 (11.88) 

SE ic (107.79) 

22.58 (8.88) 

8.2 3223) 

5.858 (2.7) 

20 (62.0) 
9.31 grams 
9.10 grams 

539.06 (23.25) 

27.18 (19. 790) 

31 - ۸6 (12.599 

22.56 (8.88) 

Tad (2.87) 

4 813 (1:895) 

16 (61.0) 


0.18 grams 
0.04 grams 





ሙ=‏ ل 


eo e 9 


Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Heignt above bed surface: cm (in) 
Static Bed Height: cm (in) 


“1: en ጠን (in Hg) 
cP: cm water (in water) 
Temp: C (F) 
SAMP # Weight of Sample 
| (grams) 
196 0. 55 
و‎ 17 0.97 
188 9.97 
8 9.94 
200 9.06 
201 9.93 
202 ፅ. 94 
203 9.07 
294 0 . 5 AVG 
205 9.08 STN DV 
Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 


Static Bed Height: cm (in) 


Pl: cm Hg (in Hg) 
GP: cm water (in water) 
Temp: ር (ዮኮ) 
SAMP + Weight of Sample 
(grams) 
206 9.91 
207 0.02 
208 9.91 
209 9.91 
219 9.91 
21 2.01 
217 0.1 
213 2.02 
214 9.01 AVG 
215 9.91 STN DV 


ہے سا لمهم - 


I uu H H MH M 


ዘ ዘ ዘ ዘ UW Hog 


53 . 5 
26.67 
32.33 
22.56 
5 

32236 
15 


59 . 5 
25.94 
35.02 
22. 56 
5.3 
L613 
16 


0 . 1 


23 
)19.50( 
اق 
)8.88( 
)2.55( 
)1.274( 
)61.9( 


grams 
grams 


(255789 
(19-25) 
(13.90) 
(8.88) 
(2.48) 
(1.0.7) 
(61.0) 


grams 


0.004 grams 





Height of 
Fressure 
Tap above 
Dretr iba POr 
(inches) 


( 
ቱ (4 ዮን ፦ SO 002 آال) (0 لد‎ & UN مہ‎ 
o) e e e e m go gm o o g o 0) 


pa pa fr fa pa 


Height of 
Pressure 
Tap above 
Distributor 
(inches) 


0(7 0( 0( 0( 0( 0( 0( 0( 0( )0( )( )0( )( )0( 


pa 
ل- 00( تی‎ MN + Q N r 


مب 
he‏ 


ኑ pa pb 
£o n 


TRRPSHEPROHT- 9412.75 in 


Pressure 
for 
Samples 
256-65 
(cm ዘ20) 


- 
48.5 
40.8 


TRAP HEIGHT = 14.25 in 


Pressure 


for 


Samples 
216-225 
(cm H20) 


52 - 2 
48.7 
40.9 
34.7 
208 
29.6 
14.9 


C O ہو‎ cC» C r9 09 
OOO OO O O) 


Pressurs 
for 
Samples 
266-275 
(cm H20) 


54.8 
49" 9 
42.1 
SS. 
23-7 
23.8 
1732 


pe 
ሠ 
کہ‎ 


Sooo 8S + 
Soo ain ۸ 


Pressure 


for 


Samples 
226-235 
(cm H20) 


534 
43. 
42. 
36. 
73. 
23 
13 
19. 


O جك‎ ጩ ጩ (1-1 =1 ሀ ጩ (በበ ھك یہ ری ٭‎ 


O O YO سا ته تی‎ 


Pressure 
for 
Samples 
2765-285 
(cm H20) 


^g! ۱ 


6 G © هت‎ (D Q) ui e ፦፡ (O ¬ + ¬+ N 


r^ rn n ol یں‎ 
بي جم به هم‎ ፡“።ር0ህ 8۵ GQ PO O |+ ት* y 


. 4 . 


Pressure 


for 
Samples 


£35245 


(cm H20) 


25.4 
51.8 
43.9 
ا 
31-7 
25-7 
18:9 
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b 
ጩጩ N 
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Pressure 
for 
Samples 
205-295 
(cm H20) 
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6 هت‎ © -J) -1 ui ^» 01 ی بت‎ 5 o n 0 


Pressure 
for 
Samples 
24969-275 
(cm H20) 


58.0 
S227 
433 
49. 
34. 
2 
21 
14. 


0 إن ہک تيت لت መሠ‏ 
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irap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
Static Bed Height: cm (in) 
PI: cm Hg (in Hg) 


GP: cm water 


Temp: 


GeCF ) 


SAMP + 


246 
247 
248 
243 
250 
281 
292 
233 
254 
cus 


(in water) 


Weight of Sample 
(grams) 


3: 
5 


34 
36 


4.64 
S. dii 
4.17 
4555 
4.55 
4.33 
8.44 


E 


20 


AVS 
STN OV 


+ value not used in computing 


Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 

Trap Height above bed surface: 
Static Bed Height: 
PIE 65۳5 ) ۱ ( 


dP: cm water 


Temp: 


Eur) 


SAMP # 


295 
231 
238 
238 
240 
241 
242 
243 
244 
243 


cm (in) 


(in water) 


cm (in) 


Weight of Sample 
(grams) 


|. 
1 
دا‎ 
1 
BS 
205 
.84 
. 89 
297 
258 


سو وم صو مسر عم رم 


S 
88 
92 
80 


AVG 
STN DV 


198 


35.20 
28.58 
7262 

22 56 
8.2 

6.398 
16.5 


(14.25) 
MIS) 
( 3.00) 
(8.88) 
(3223) 
(2278) 
(62.0) 


4.48 grams 
0.72 grams 


average. 


36.20 
2130 
8.90 
2255 
9 

4.768 
1! 5 


(14.25) 
(10.75) 
( 3.50) 
(8.88) 
(2.87) 
Clos) 
(52.0) 


1.88 grams 


0.22 


grams 


 . - 


P, 
TIMLULA 


5۳ 


Y 





Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
Static Bed Height: cm (in) 


Pi: € Hg (in He) 
GP: cm water (in water) 
Temp: ር (የዮ) 
SAMP E Weight of Sample 
- im (grams) 
225 0.93 
227 0.68 
228 0.61 
229 9.56 
230 0.81 
231 0.52 
232 0.38 
23 2.5 
224 0.75 AVG 
235 0.44 STN DV 


Trap Height above Distributor: 
Bed Height above Distributor: 
Trap Height above bed surface: 
Static Bed Height: cm (in) 


cm (in) 
cm (in) 
cm (in) 


Pl: cm Hg (in Hg) 
GP: cm water (in water) 
Temp: C (F) 
SAMP # Weight of Sample 
(grams) 
216 0.32 
217 0.35 
218 9.36 
209 0.30 
220 0.32 
221 2. 25 
222 90.24 
ረጋ 90.30 
224 0 . 4 AVG 
225 0.30 STN DV 


y 
IÓ 
(ይ 


ዘ og 


25.26 
26.87 
2.53 
22.56 
6.6 
3,218 


18.5 


nn uU ክህ ዘ Ul 


14225) 
(10.50) 
(3:73) 
(8.88) 
(2.560) 
(12267) 
(62.0) 


0.52 grams 
9.17 grams 


356.20 
25.04 
10.18 
22.56 
3 
2.524 
15:5 


ዘ Uu ዘ ከ [[ "n Ht 


= 9.32 


)14.25( 
(10723) 
(4.00) 
(8.88) 
(2.48) 
(1.0933) 
(62.0) 


grams 
0.04 grams 





32397 0 
28.58 (11.25) 
A281 (1.50) 
22.56. (9.88) 
8.1 (319) 
6:135 2 (2775) 


Trap Height above Distributor: cm (in) 
Bed Heıgnt above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
Static Bed Height: cm (in) 

Pie Gh He. Cin Hg) 

GP: cm water (in water) 


Temp: C (F) 17 (83.8) 
SAMP 5 Weight of Sample 
e E (grams) 
286 = መ 
287 Tas 
288 7.46 
289 8.05 
299 62 
Z34 1-26 
292 ርን 
228 7.855 
294 7.04 AVG = 7.22 grams 
ey 295 == 648 STN DV = 8.51 grams 


3 ۹پ‎ 
27.30 (10.75) 
5.03 (2.90) 
22.56; (8.88) 
p (2.87) 
4.836 (1.304) 
17 (53.0) 


Trap Height above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 
Static Bed Height: cm (in) 

Bi CE ES (in Ag) 

dP: cn water (in water) 

Temp: C (F) 


' SaMP # Weight of Sample 
(grams) 
LG 2.44 
2 7 3-7-7 
278 3.58 
23 4.04 
280 4,56 
له‎ Sat 
292 4 کی‎ 
299 5.53 
234 S375 AVG = 4.35 grams 
285 5.49 STN DV = 1.05 grams 








Trap Height above Distributor: cm (in) 
Bed Heignt above Distributor: cm (in) 


Trap Heignt above bed surface: cm (in) 
Static Bed Height: cm (in) 
Pi: cm Hg (in Hg) 
GP: cm water (in water) 
Temp: C (F) 
p SAMP $ Weight of Sample 
(grams) 

266 9.74 

267 1.80 

268 1. 7 

263 fal 

67 . اا 270 

271 15 

21 1.8 

72 1 3 

274 1.60 AVG 

25 1.6 STN DV 
Trap Heignt above Distributor: cm (in) 
Bed Height above Distributor: cm (in) 
Trap Height above bed surface: cm (in) 


Static Bed Height: cm (in) 


Pl: cm Hg (in Hg) 
dP: cm water (in water) 
Temp: ር (ዮ) 
SAMP # Weight of Sampie 
(grams) 
256 0.59 
257 0.76 
258 0.67 
£58 0 5 
269 8.54 
251 0.52 
262 8.47 
263 0.65 
264 90.51 AVG 
265 0.87 STN DV 


rJ 
© 
p 


23209 
29 - 
Bee 
222.96 
6.6 
3 5 
15.3 


1ھ 


)1 m) 
) 
2,28) 
(8.88) 
(2.60) 
(1.287) 
(62.0) 


grams 


0.37 grams 


32.33 
26.04 
5 
22. 56 
53 
2.684 
18.5 


ere i=), 
019225) 
(2,50) 
(8.88) 
(2.48) 
(1.025) 
(62.0) 


0.63 grams 
0.13 grams 





J‏ لا یح حم مم 


This Appendix contains the output from the image analyzer. 


ine output is arranged in the following order. 


Bed distribution 

4 cm freeboard height 
8 cm freeboard heignt 
tZ cm freeboard height 
18 cm freeboard height 
22 cm freeboard height 


31 cm freeboard height 


Each of the distributions at a given bed height applies only 


to the Uo/Umf = 3.81 condition. For each of the seven image 


analyzer outputs listed, the following three formats are used. 


انا 
s‏ 


Histogram of absolute particle frequency 


vs particle diameter. 


Cumulative percentage plot of particle 


sizas distr:oxttion. 


Table listing of the above data. 


202 





The dashed line on both the bar graphs and on the cumulative 
Percentage plots represent the Geussian distributions which fit 
the given set of data. This Gaussian distribution should be used 
only as a rough estimate of the data because the analysis included 
“particles” less than 70 microns. By viewing a blank slide with 
only the tape aeplied, these “partiicles" were confirmed to be 
busoles end Girt entrapped in the adhesive on the tape used to 


hold the sample particles in place. 
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Height Distribution 
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APPENDIX K 


This Appendix contains the photographs of the oscilloscope 
traces obtained while sampling. Each of the pictures is labeled 


with the sample numoder for which it represents. 
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This Appendix contains a complete listing of the particle 
trajectory computer mode! used in this research. In addition, an 


in depth program description and flow diagram is included. 


نم 
ol‏ 
ርበ‏ 





Model Computer Program 


This section describes the computer program listed in 


Appendix KL. 


To assist in understanding the program logic, the 


program itself contains numerous comment statements. The program 


consists of a main program and seven (7) subrcutines. The main 


Program controls the input of parameters, trajectory calculations, 


and output selections. The subroutines control the actual data 


output in either table form or graphics. The following is an in- 


depth description of the program: 


Lines 


58-158 


122-610 


P. E: Description 


1 


This section is used to explicitly define 
the major variables used within the program. 
The large arrays are defined in common block 
form to save memory.  Constants used in the 
program are also defined. 


Tnis section is used to input variable data 
to the program. Three options exist for the 
bed particle distribution input: 


1) The default condition sets the quantity 
of each particle to unity. This option is 
Generally used when a height determination 
15 required or being sought after. 


2) The particle number option allows the 
entry of bed distribution by the number of 
particles present in each diameter range. 
This is used when cata from the image 
analyzer was being used. 


3) The last option is used when the particle 
size distribution is determined from a sieve 
analysis and the data is measured in grams 
mass. The program will then determine the 


256 





558-530 
700-710 


720-800 


819-1659 


829-849 


853-1340 


869-3949 


1010-1060 


1100-1170 


1180-1220 


1289-1320. 


1358-10 


1380-1430 


1440-1500 


particle number density based on the 
assumption that the particls are spherical. 


These lines convert the velocity input 
values from Ft/s to m/s. 


Determine the sloge of the triangular jet 
using the input amplitude and duration. 


Initialize the date arrays to zero (0). 


This section calculates tne trajectory for 
each diameter particle. This WHILE condition 
contains the following 8 subsections. 


Initialize the height of the particle ta 
pass the first WHILE statement. Calculate 
the particle diameter to be used. 


Calculate the particles trajectory 
parameters while the particle is above the 
bed surface. This DO loop contains the 
following 6& subsections. 


If the elapsed time since the particle left 
the bed is less then the jet ducation time, 
the add the jet velocity to Uo. 


Calculates the relative velocity of the air 
with respect to the particle. Determine the 
sign of the drag force. 


Calculate Reynold’s number. 


Calculate the drag on the particle and then 
determine the particles neu velocity and 
position. 


If the particles velocity is positive, 
record the new max height and elapsed time. 


Add 1 count to the probability array in the 
storage position representing the particles 
height/2. By using the DIV statement, a 
height window (dH) or Z cm is created. 


Change all heights from Ft/s to m/s 


Multiply probability distribution by bed 
distribution weighting factor. 


Determine the maximum value and normalize 
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1510-1569 


1578-1660 


2268-2520 


2560-39160 


2540-3129 


3130-3240 


3250-3830 


3840-3928 


3380-4130 


4130-4289 


4330-4468 


4510-4648 


the data such that the maximum is 100. 


Calculate the narticle entrainment by 
summing the total volume of each particle 
distribution at a given height. 


Normalize the entrainment data to 19% and 
calculate the ln value if not equal to zero. 


Display outzut menu. 

Controls the selected output. 

End of MAIN program 

Sub Disolay_Data is used to list input 
perameters used and the maximum height 
attained by each diameter particle. 

Sub Display_graphics controls the processing 
and output of graphic information. This 


subprogram conteins the following 4 sections. 


Performs scaling for determining graphic 
dimensions and limits. 


Creates graphic display and labels the X and 
Y axis. 


1 


Controls the plotting of selected graphic 
output. 


Prints hard cooy. 


Function routine for determining scaling 
factor for graphic dispiay. 


Finds the limit of data within an INTEGER 
array to limit the X axis on a plot to the 


range in which the Y values are non-zero. 


Determines the maximum value stored in an 


INTEGER array. 
Same as 4330-4460 except for REAL array. 
Same as 41886-4280 except for REAL array. 


END of PROGRAM 


IN 
ما‎ 
0 





12 


29 | 

39 | 

43 | 

20 OPTION 3asz 1 

22 INTE Answer ‚Soec:fy,Point ‚Vieu,Bed_nın, 1 

70 ውሪ ፣ ንር لی‎ +, Part heısnt ,Slooe, Jat.t,Jt 

29 Can 1INTeser Disiribution(39,90),Distr:ib_dens1iy(33,80),REAL Height(2,50),3 
sd gpar:icisQs0).tat72inment(2,80) 


MAIN Program....... 


32 DIM Somes(21,Clears(7] 

190 Clsar SC 223( 123 JCA2RS(7<) ! CLEAR SCR kay 

119 Homes*CHRS( 255 )&CARS( 8+) ! HOME key 

129 Uiscos:ty xine1.i188c-z I air, m*m/s 

139 Densıty_aır=1.201 I kgim+m®m 

140 Densiiy_part=7.362+3 ! kG/memem 

1S9 Graviiysz.3085z2 | m/s*s 

159 ! 

179 | Incut var:iaoles 

¿32 l 

190 SQINTZ2 TS 1 ١ 01۴3٣1 to CAT 

ae INPUT “Enter mean air velocity Ue (Ft/s): *,Ua 

200 ۲۱۳ ال‎ enter amplituce of jet velocity (Ft/s): ",Ui 

EZ INPUT “enter initial particle 0 Uso (Ft/s): *,Usol 

230 101۳۱ 23۵ duration of jet (s): *,Jet_: 

256 - - INPUT “Enter time increment for Se aa (32: ",Del t 

252 INPUT "Do you desirs to input tne Bed particle distribution (1 yes,% no): 
,Ansuer 

260 TF Answer=@ THEN --- | Reg mist cant wastes f Bed dist not wanted 
270 FOR I=1 TO 9۵ | Set bed dist to [1 
559 Sed 3articie(I)s!. 

E NEXT I 

209 EE 

310 IF Answer=1 THEN ! Bed dist wanted 
322 INPUT “Inout will be; weight in grams (1) or numoer (2): “,Scecify 
330 IF Soec:’v=1l OR Soecify=2 THEN 

342 PRINT "Enter data for sach diameter" ! Enter gram or + 
3350 Min_b=d=1 

26% POR i=! TO ۵ | Sach diameter 

372 PRINT سن"‎ .: ۰10+ ۲۵ ۹ 

>89 INPUT Bed 3article(CI) 

220 IF Sed_particla(I )<Bed_particle(Min_bed) THEN Min_bea=[ 

402 NEXT I 

510 ፲ኛ Soecify=! THEN | 1? ueignt entry 
120 FOR T=! TO 0 | Calculate + 

4 0 Sed _sarticle(!)*6*8ed_particle(i)/(Pl+(1T*10+70)*5*Density_nart ) 
552 IF Sec 6ه 1(8 )هد 601:: م23‎ nariiclae(Min Beg) THEN Min bed-I 

452 NEXT [ 

+62 ENO IF 

273 EST 

£63 eons 220 

0 ND IF 

229 Sa 

S10 SO10 250 

2-9 END IF 

< 20 END IF 

549 Max_bea=! 

==22 FOR ۲۶۱ TO z2 | Sind largest + 
220 TF 3eg 2ar:icle([IYX8eg nar::cle(Ma« 3egd) THEN Mas, 22271 
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سڈ ۔ ۱١‏ و داء : ە دلا 
2*28 ۰ ۱ اعد دلا 
21 27ل YES jet ii THEN‏ 


37 NEXT I 

36 Max_value=!90/220c_cart:icla(Mas _ hed) ! determine scaier 
33 FOR fet TO 30 Races zata 2 ۶ 0 
29 Sec_yartıcial.!=Beo_ Sarııziell)eMax_value 

31 NEXS I 

5ረ ۱ 

E3 ! Comcute data 

2= ۱ 

23 U12Uo 

1 2۱ علا 29 

27 Us tijoe. 3928 | Change to m/s 

5 

7 


Jo C2 44) በ) 


Oa Ya Qa نیم‎ ር3 ሮጋ ሮ3 ወ C9 C) رم‎ ዑሮ3ጋ ዕህ ሮጋ ሮጋ C0 C) C) م‎ C) CD C) C) C9 (0 C9 C2 9 C2 


Siwger2.J pp/det t ! Triangle z2t pulse 
73 FOR i21 10 zà 
13 He:ant(1,1)28. ! Zaro height array 
75 حا‎ > ) 2. ۶ ۰ I Zara numser array 
7S SOR J-t TO 39 
73 DistasDetronl i, J 8 | Zero disi erray 
7a کو دی روج ےد ہد‎ Üi)=0. | Zerg antrain array 
78 =aAtrarnment(2,J)=d. 
7 NEXT J 
39 NEXT i 
31 -GR i=! TO 0 | Cer each perticl= 
9 FPart_nergnt=. 38080! ! set Tor WHILE stats. 
23 UcosUcol | Initial V constant 
3£ Diamere”_aria(i-7)e1.03-= ! en Dia um 
ez UHILI Part heigni»Q 
3s مدع زک چاو‎ 
$7 1i*e nousMe:ognt:(2,1)e0281 5 !ا‎ lims of nern 2116+ 
38 T= Tıme_now(Jet_: THEN IT Ber: ڑج‎ 7377 
225 2 ሆረ لے ال ,2,۱ ) کال‎ ۰+ ۷ | 1st hal? of jet 
2909 Uo*Uo-Sioge*rde:gn:(2,:)*50e1l i 
50 == : 
328 ۳02و : 1-51006۰۳2 ز لا+ول۱<ع)‎ ) 2 , (۰۱ 1 _ + ! 2nd half of jet 
330 END ir 
2:9 zNO I? 
229 Aelative velsUo-Uoco ! Rel U seen Dy part 
2252 3ign=1 | Set «iom positive 
272 se Relative_vei<d THEN I is V negative? 
coe ۳ [1 ! Sei sign negative 
sog WE :፦ 
1223 ۱ Calcuiate ٥ 
1419 Aeynolcs no*8eiaiive vel*Üiamseierc pari/Uiscosiiy kin 
1920 Revnclas nosaB8S(Reynolas no) 
1929 ir Reynolds nosQ. THEN 
1242 eacaleratione-sravity 
1959 GOTO 1159 
1953 ENO iF 
1979 | Caiculate drag 
1089 il Usas data correlation for drag cos??:2:2ni gooc for > 
1498 ۱ 
1130 Drag coef f=24/Raynalss_notd/( 1+SOR(Reynolds_nc))+.4 
1 1 ÜragezigneUm-ag zoeíf*Density ai7*Ceiative vel*Raelative vele7:*0iameier 
.2erl*Ülameiec pac:/&8 
1119 ١ veiculats acceleration 
1138 Acceberas:ometrac/ (Sens: ty_sart +? i eh: amece- sart*s/c i-sravity 
1142 ١ Caiculats gal velocity 
1199 UelociiysUZ2oT^ccelerat:io0neÜei - 
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Po of h ehib 14 444 


C) 1460 ዐን -1 ብ (በኦ OI fie O in O1 = 


m Im سر سو‎ œe œ pa la ya po po o yoo 


ዐ) -1 ዐ) (1 LI 63 c‏ 19( دک 


سر میم quo lt.‏ مسج ye be‏ سم مہم 


=p =p =) >=) ۰۱ لل‎ ۱١١ (010 (11 ۱۱۱ 09 0) on O) الغ‎ O MN AN ۲۱ ۱۶۱ 4/۱ uU (Y UV 0 (Y 9 & ها حغ‎ 


۱ ۱۱۱ (1.16ሁ፦ 


(በ fe ا٠‎ ۰1 e“ لځ‎ ved ae 


سر fr y. (1 p p2‏ مام سیا 


l Calculata new position 
Part_he:gnt=Part_hergntt(Velocıty+Uno )+-Del_t/?2 


IF Upo>@ THEN | Part still rising? 
Hezgnt(1,1)=Part_hezgnt I Save position 
Hergnt(2,1I)=He:zgnt(2,3)+1 ! Inc t for max rise 

SNO IF 

Uoo=Velocıty ! Set new Y for next inc 

Pornt=(100*Part_nezgnt) DIV 2 I 2 cm wıde storage bins 

I7 Paint» 79 THEN I Set default for fatal 
Poınt=?2 

ZNO ۳ 

IF Poin:i<=0 TEEN ! Set default for fatal 
2۵0۵9 

ZNO IF 

PRINT T,Point | Indicate como working 

! Save + times part in neignt bin 


Disiribution(l,Po:znmt+l)=Distributionti,Point+l)+1 
=NO WHILE 
NEXT I 
FOR اع‎ 10 Se 


ሣመ ሮጠፒ ( ፲ , ፡፲ ۱۶۳۵106) 1 , 1 ۰ 089 ! Change to cn 
NEZI 1 
FOR [#1 10 ۵ 
Ualuszdec 2naciicletI) l! Part siza weight factor îron hed dist 
FOR J=! TO 30 
Distirib_aens:ity(i,J)=0istributionti,J)*Yaelue ! Weight dist values 
NEXT J 
NEXT I 
Mex_dist="NMax_int(Distrib_density(+),50,80) I Find max value 
Factoc-100./Max dist I’ Scale for 100 max 


FOR 13۷ ۱0 0 
FOR J=1 TO 99 


01 52:5 6۳7511۷ ,[ (2: 5۳:0 32051۷) ؟‎ , [ (۰۴ 6+07 ! Scale values 
NEXT J 
NEXT I 
FOR I=! YO SQ I Mass densıty/unit area at neignt aoove bed 
VolumezPT+((1-7)/1000)73/6 I Volume Cu cm 


FOR J=1 TO 0 
entrainment(l1,J)=Entrainment(1,J)+Volume-»Distrib_density(l,J) 
NEXT J 


NEXT I 
۳13: ۶7۲۳ 5 
FOR I-1 TO 30 ! Find maximum 


IE Satraınment(1,I)>Entraınment(l,Mar_entraın) THEN Max_entraln=l 
NEXT I 
Pactor=!00/Entrainment( 1 ,Max_2nterain) ! Normalize to 0 
FOR ፲=! ፲0 398 

Entraınment(L,!)=Äntraınment(|L,!)>Factor 

ir entrainment(!,1)«-Q0. THEN 4 

znrraınment(],!)=L0Giintrainment(!,2)) 


NEXT I 

! 

e gnt^ol 52-8 

| 

29 0/0۳۶4 IS i ا‎ 01,127 to CRT 

USING "5/"‏ یی اس 

PRINT "1) Displey hezgnt vs diameter data” 

POINT "2) Dısolay ner:gnt vs diameter grann” 

gern, “3? Bol, ensity ysihe:gnt 42s function af dia graon” 

SCINT ^-^ Pisolay tensity vs diameter as a function of heignt grcaon 
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NJ 
O 
ወጠ 
O 


2070 
2080 
2928 
2190 


21'0 


1) 
„> 
1) 
ርጋ 


عدو 
C4‏ 
ጩ‏ 


P) يخم‎ (3 1ህ የ) mw ۰ par 
C (0 OM =) O ÛU) te 
C يج هت ىم ےت‎ © 


ld ja هو مر صمو‎ jos y 


EN 


int 57 Disolay same as 3 but with bed density’ 
INT ^3) O:solay same as 4 out with oed density” 
INT ^7) Disolay density vs diameter of bed mass” 
INT ^3) Disolay entrazmment density above bed” 
INT *9) Disolay Ln antraınment densıty anove bed“ 
INT 9 EXIT PROGRAM? 


PUT ^£nterc numbec of desirsd display: ^,Ansuec 

LEITT Answer 

CASE =! 

CALL Oisolav_data(UVi,Uci,Uj,Jet_t,Dei_t,Home5,Clears) 
= وې‎ 


= 

CALL Disclay _3zjraoní(CizarS,Aomea3, Answer, Vizw) 
s 

T 


INPUT “Enter و-< رو عاح: مود‎ to be vieued (80-70 um)(0 foc all): ",Uieu 
IF Viaw= THEN 0 

IF View<83 OR View>S7Q@ TREN 0٥ 

Uiazus(Uisu DIV 10)-7 


CALL Disclay_graoníClears,Home3,Aánsuer,Vicu) 


CASE = 
INPUT “Enter dasıreo heignt above bed surface (93-153 cm): "Viau 
Tr View<S OR Vizw>153 TEEN 1953 
Uieus(Uisu OIU 22-41 
CALL DOisolay_graoní(Cisars,omess,ánsuer,Viau) 
۳: >ہےہ‎ 
INPUT “Enter particle sıze to be viewed (39-3798 um)(®d for all): " ‚View 


IF Viaw=8 THEN 2339 
TF Viaw<82 OR Vieuw>S79 THEN 2000 
Vieuw=(Vizw OLY ۱۵ (-7 
CALL Oisaelay_grapníClear3,Home3,Ánsuer View) 
CASE -5 ^ 
INPUT “Enter desired heignt above Bed sucface (Q-158 cm): "Visu 
IF ViewxQ OR Viaw>1S8 THEN 2269 
View=(View OLV 2)+1 
CALL Display_graon(Clears, Somes ,aAnswer ,Viau ) 


Case -7 

CALL Or.solay_sgraan( Clears ,HomeS answer, View) 
60555 ٦ 

CALL Display_graon(CiearS,Homes, Answer ‚View) 
CASE =3 


CALL Display_z3raoní(CizarS,Home35,Áánsuwer ,Vieu) 
65 10 

STOR 
CaS. SSE 


GOTO 18292 


Dose مات‎ 


2210 GOTO 0۵ 


¿220 END 

2230 |! 

ci |! 

2 20 |l 

2-20 ' Sus usad to lis? inout ceramatars and mar ۷م‎ 55 at time + ser diameter 
75 ۲ ZU: 02 23139. atal Ui Jol, Un, Jai 2 Mei: somes. clears) 

«ዳሸን... ርን 1”:0ህ፣ወሰ(ቅ/,:5፡:::ሮ. :መጣ5:2ሃ(”።7ያ 5۹ He:ignt(*),S3ed narcicle 
)د ود و ونت ږل‎ 

۲۳ 2۵ 0:281 ت1 5 اا‎ I Home and Ciar screen 
2:00 OUTPUT 2:Clears: 

2 7 ! Prınt output 

2 ¥ PQINT ^" cm/s" 

2230 POINT “Mean 223 Nalacıty= س۰۰۶٥‎ 

223-09 FEIN, “initzai Part:it2 VYelgeıtv= ولا‎ ۱۰ 2 . 9 





2359 PRINT “Pesak Jet Valoc:ty= " : ز لا‎ ۰ 30 . 8 


2360 PRINT “Gas Jet Ouration= e t 

2378 PRINT 

2380 PRINT ^ Diametar”,” Heignt^,^ Tine", Diameter”;” Height^ 
Time” 

2290 PRINT ° um,” ar ^ Secs vs um^,^ cmo T seconds 
2490 PRINT 

¿2210 FOR I=! TO 58 STEP 2 

2120 EM USING 20=%,20,5X,30.09,3X,0.000,5% 3°: 1#7 3)» ,Heschtt L, I), Heightt2 
some! t.C0I98)0€910,He1g9gnr(l,I91),Heignt(2,I1*1)«081. د‎ 

2239 NEXT I 

2442 PRINTER Se! 

2553 (۲ mount hara copy? (l1)= yes, (9)= no: *,Answer 

2463 f= Answer=! THEN 

2470 PRINTER iS 701 

2480 59ፐ0 22260 

جع ے2 2420 

2590 TF Answer<>@ THEN 2452 

25109 =NO IF 

2520 SUEEND 

25360 ! 

2328 |! 

¿338 | 

2550 ! Sub used to control graonics outout of data 


حم هم ممم 


¿579 SUB Disolaey graon(Cleacs,Homes,INTEGER Data_set,Vizw) 

2588 COM INYESER Disiribution( +), Distrio_density(+),REAL Hergnt(*),8ed_ sarticle 
(+) Entraimment( (ء‎ 

2530 - REAL Xmax, Ymax ,Xtick Ytick Xmin,Ymın 
¿690 OUTPUT 2;Homes; 

2619 OUTPUT 2:Clears: 


iU 


2529 SINIT ! Initialize graohics 
2539 GRAPHICS ON 

2590 SELECT Data_3et ! Scale plot routines 
2539 CASE =2 

¿500 Xmax=500. 

2570 Ymex=FNMax_rzaal(Heignt(+),1,50) 

2568 Xt1ck==4. 

50 YtickzrNScale( Ymnax) 

2109 Case 5 

2710 Ymax="NMax_intí(Distributiont+),24,80) 

2720 Xmax=PNData_!limit(Disiribution(»),54,30) 
2720 tıck=rNScala(lXmax ) 

2740 Ytick=rNScala(\Ymax ) 

2752 CASE =4 

2753 Xmax=s0@. 

2779 Ymex=FNMax_int(Distrısution(«), 59,39) 

2780 ጀኒ :ርኗ==26. 

> 20 Yııck=rNScale(Ymax ) 

2800 Case =5 

2919 ۸62 ۶۶۲۱02 ۱۵ limit(Qistc:b dens::y(C*2,229,30) 
99 :፡ ۷2 2۶ 0 ۰ 

( دب ۸ )2 ددع ረረ.‏ 2230 

ز رد )دز دم حول رن د۷ sa‏ 

2899 CASE == 

2960 Xmax==30. 

>= mg Ynarz-NMax int(Oistp7:b zensitiy(*)2,:22,30) 
¿380 Xtick=850. 

2282 Ytick=FNScalal(Ymayx ) 

> 2930 CASE 7 


Dub 
fs 
Il 
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2566 


Xmax2500. 

Ymax=!00. 
tick=sd. 

Ytıck=19 

CASE =3 
Xmax*-NOata limitr(Entcainment(*)2,1,30) 
Ynax*120. 

Xtıck=FNScal=(Xmax ) 
ECE. 

CASE 3 
Xnax="NData_limıtr(Entrainmentt +»), 1,39) 
Ymax==. 

Ymin=-3. 
Xtick="NScal=z(Xmax ) 
Ktick=,5 

የ ከገጹ ፡ዓ።ሙ። = 

SOTO 3112 


Ymins-Yiick 
Xmins-2.*Xtick 
X 3ጂ፳25=5. 


UINDCU Xmin,i.l1*Xmax,Ymin,l1.l1e*Ymax I Set graohics scale 
AKER uck.Ytick.X axis.Y axis ! Set graan axis 
LORG S ' Lagel X-axis 


FOR [=X_axis TO Xmax STEP 2 
MOVE I,Y_axis 
LABEL I 
NEXT I ٢ 5 
LORG 8 ¦! Label Y-axis 
FOR I=Y_axis TO Ymax STE? Ytick 
HOUE X axis,I 
۳ 
NEXT I 
Sebel: Data. sad ! Plot data 
CASE 2 
MOVE 80 ۳۶: 33۶41, ( 
FOR f=2 10 90 
ORAYW (TI+7)+19 ,Hezgnt(1,1) 
NEXT I 
CASE 5 
Segin=View 
Pinisn=View 
IF View<2 THEN ۰ ٣۴ 9 
Begın=! 
Finisn=S0 
>0 ۴ 
FOR [=Segin TO Finish 
MOVE sO dat” ibut:ont 1,1) 
FoR J=! TO ۷۸ 2 
ጋመ ۰۲5 buU onli.) 


zx 
-1 
۰ 


ጋ 


> 
LI 


۱ - ۷ ور 
i»‏ 


(=s; 


p 
5 iN a 
€) )3 It 
= 
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A) 


(3 


erica. Una) 


AY > 
~) 


«10, istribution(I,U:eu) 
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i 
۱ 
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C) C) C) ين دخ‎ 2 C20 ኔኔ ኔር 9 C20 C. C ی مت‎ ሠጩሠ 


al => ل- 1= 1“- لل ام اه‎ OI 0) 0( OI dai OI dg) a d) aiu 


۷ بح 0 O =} Gi (ዘብ ሁ 04 PL‏ قي O‏ ہام لام یں UB (N f‏ لد 


۱ 604 O} OF GF GE OF OF Gi bbb GE GE GE Ul OG Ot OF OF LE © 


د 


وہ 
ጩ‏ 
ጩ‏ 


pr 


tinishzUieu 
IF Uieu«Q THEN 
2=2 : በ=! 
Finisn==U 
END IF 
FOR I=3esın TO Finisn 
MOVE 3, Distr:ib_densityíl,1) 
FOR Jz=! TO Xmax/2 
DRAaU 2» J Distria_aensitytI,J) 
NEXT J 
NEXT I 
CASE =3 
HOUESSOSDISICSE density(! ,Vieu) 
FOR I=2 TÛ 3@ 


! 21862 7080 


DRAU (I+7)+19 ,Distrib_densityti, View) 


NEXT I 
CASZ =7 
meve 38 9 6۵ :م2 مر‎ ۶) 1 ( 
FOR 12 10 ۵ 
ORAW 1۲۰۱۵ , 322 +حو مر‎ 61) ( 
NEXT I 
CASE =3 
MOVE 0۵ , 20:۵۵ 
FOR I=2 TO 0۵ 
DRAW [+2 ,Entirainmentí(1,I) 


NEXT I 

Case =9 
MOVE 2,Zntrainment(Z,1) 
FOR I-2 TQ 80 


IF Entrainment(2,1)=9. THEN 3849 
DRAW [e2 Entrarnment(Z2,1) 
NEXT I1 
5۳ 2 
INPUT 32۳8 graon (1) ves, 
TF Answer=! THEN 
DUMP DEVICE IS 
DUMP GRAPHICS 
SCLEAR 


(4) na: 


401 


IF Ansuerc«»0 THEN 3840 
=NO IF 
50-58 
SUBENO 
! 
١ 
| 
| Function determines qgragnics axis 
DE? 7NScals(Ymaz) 
IF Ymax»100 THEN 


(፤ 


THEN‏ 22 روم۲ 2010 2۵0 لح ده 


* Answer 


I! Cumo to printer 


scaleing 
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END IF 
RETURN Tick 
FNENO 
| 
l 
۱ 
۲ Function detesfmimes mex extent of date for X-axis limit 
=> FNData_lımit(iNTESER Data_l(*),Row_mar,Col_nmax) 
521=1 
REPSAT 

Sum=3. 

FÜR i=1 TO Row_nax 

Sum=Sumrdata_1(1,Col) 

NEXT I 

Col=col>1 
UNTIL Sum=Q 0% Col=C91_max-1 
aSTURN 2۵ 
N 


[ 

[ 

1 

t Function determines max value in integer array 

05= FNMax_intt INTZSER Data_1(*},Rew_max,Col _max ) 
ጠመ! 


OR I=1 TO Row_max 
FOR J=1 TO Col_max 
IF Bata 1(I,J)2Q0ata !(R max,C, max) THEN 
፳ maxzI ١ 
C_mex=J 
END IF 
NEXT I 
NEXT I 
Max=Data_l(R_nax,C_max ) 
RETURN Max 
FNENO 
١ 
! 
۱ 
I Funczion determınes maxımum value in real array 
QE? FNMax rsal(Qata 1(*),INTESER Rou max,Col max) 
R_nax=l 
C_max=1 
FOR I=! TO Row_mex 
FOR J=1 TO Col_max 
IF Oata_1(1.J)>O0ata_!(R_max,C_max) THEN 
R_max=! 
C_mex=J 
END IF 


(AMpex, Lomax ) 


11 

| 
E 
aA) 
22 


5 twee 


NENO 


unction determines max 2xitant of data for X-axis limit 
-*NData limitr(Qata 1(*),830u, Col max) 
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IF Data i(Rou,Col)4s0 OR CalzCol max THEN 
Limit=tole2 
RETURN Limit 

ENO ir 

Col=foi+1 

GOTO 5716 

FNEND 
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This Appendix contains the calibration data for the 
anemometer probe. The calibration was conducted in a smal! wind 
tunnel using a pitot tube connected to a micromanometer capable of 


measuring pressures to within 0.001 ins. of water. 


CO 


۱۰ 
Ja 





Air Velocity (ft/s) 





9 1.0 21 3.0 


Oscilloscope Trace Amplitude (Volis) 


FIS. M-I Calibration of Anemometer Probe. 
Oscilloscope Voltage vs Air Velocity. 
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